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     

Mach Wave Radiation by Mixing Layers. Part II: Analysis of the Source Field

In Part II a large-scale sound source in a time-developing planar free mixing layer is studied using an acoustic analogy approach. It is shown that only the non-compact character of the source resembles well the character of the corresponding source in a spatially developing flow. A model based on a continuous assembly of wave packets is derived and applied to direct numerical simulation results of two supersonic time-developing mixing layers undergoing transition to turbulence. The analysis predicts two distinctive dominant Mach wave sources in agreement with the direct analysis of Part~I. The first dominates during the stage of the Λ-vortex structure and the second just prior to the final breakdown to a fine-scale structure. The convective velocity of the second Mach wave source is higher than the first and thus its Mach angle of radiation is higher. The second source has a reduced strength at the higher free-stream Mach number. Directivity and frequency spectra compare well with the results of Part I, demonstrating that the assumptions inherent in the analogy are quite reasonable. Received 5 August 1997 and accepted 6 April 1998     

A Conceptual Study of Cavity Aeroacoustics Control Using Porous Media Inserts

In this study, an integrated flow simulation and aeroacoustics prediction methodology is applied to testing a sound control technique using porous inserts in an open cavity. Large eddy simulation (LES) combined with a three-dimensional Ffowcs Williams–Hawkings (FW–H) acoustic analogy is employed to predict the flow field, the acoustic sources and the sound radiation. The Darcy pressure – velocity law is applied to conceptually mimic the effect of porous media placed on the cavity floor and/or rear wall. Consequently, flow in the cavity could locally move in or out through these porous walls, depending on the local pressure differences. LES with “standard” subgrid-scale models for compressible flow is carried out to simulate the flow field covering the sound source and near fields, and the fully three-dimensional FW–H acoustic analogy is used to predict the sound field. The numerical results show that applying the conceptual porous media on cavity floor and/or rear wall could decrease the pressure fluctuations in the cavity and the sound pressure level in the far field. The amplitudes of the dominant oscillations (Rossiter modes) are suppressed and their frequencies are slightly modified. The dominant sound source is the transverse dipole term, which is significantly reduced due to the porous walls. As a result, the sound pressure in the far field is also suppressed. The preliminary study reveals that using porous-inserts is a promising technology for flow and sound radiation control.     

Noise radiated by a non-isothermal,temporal mixing layer. Part I: Direct computation and prediction using compressible DNS

Direct numerical simulations (DNSs) are performed in order to study acoustic emissions generated during the transition of isothermal and non-isothermal mixing layers. The sound from temporally evolving mixing layers is computed directly using DNS for a computational domain, which includes both aerodynamic and acoustic fields. Good precision of the computed acoustic field is ensured by using a numerical code based on high-order finite difference schemes of quasi-spectral accuracy. Two- and three-dimensional simulations of mixing layers are performed for various Mach numbers and temperature ratios. For each case, the acoustic radiation of the mixing layer transition is investigated. Comparisons illustrate the importance of the combined effects of temperature and Mach number on the acoustic intensity. Qualitative agreement with existing experimental observations for hot jet flows is observed. It is also found that the appearance of three-dimensional motion leads to a substantial reduction of sound emissions. In the second part of this study, DNS data are used to perform acoustic analogy predictions. Excellent agreement between direct computations and predictions is obtained in all cases. Analysis of the source terms yields a new interpretation of temperature and Mach number effects, based on the predominance of one term over the other.     

The radiated noise from isotropic turbulence

The noise radiated from isotropic turbulence at low Mach numbers and high Reynolds numbers, as derived by Proudman (1952), was the first application of Lighthill's Theory of Aerodynamic Noise to a complete flow field. The theory presented by Proudman involves the assumption of the neglect of retarded-time differences and so replaces the second-order retarded-time and space covariance of Lighthill's stress tensor, T _ij , and in particular its second time derivative, by the equivalent simultaneous covariance. This assumption is a valid approximation in the derivation of the ² T _ij /t ² covariance at low Mach numbers, but is not justified when that covariance is reduced to the sum of products of the time derivatives of equivalent second-order velocity covariances as required when Gaussian statistics are assumed. When these assumptions are removed the changes to the analysis are substantial, but the change in the numerical result for the total acoustic power is small.This paper is based on an alternative analysis which does not neglect retarded times. It makes use of the Lighthill relationship, whereby the fourth-order T _ij retarded-time covariance is evaluated from the square of a similar second-order covariance, which is assumed known. In this derivation no statistical assumptions are involved. This result, using distributions for the second-order space-time velocity squared covariance based on the Direct Numerical Simulation (DNS) results of both Sarkar and Hussaini (1993) and Dubois (1993), is compared with a re-evaluation of Proudman's original model. These results are then compared with the sound power derived from a phenomenological model based on simple approximations to the retarded-time/space covariance of T _xx . Finally, the recent numerical solutions of Sarkar and Hussaini (1993) for the acoustic power are compared with the results obtained from the analytic solutions.     

Balanced Atmosphere–Ocean Dynamics, Generalized Lighthill Radiation, and the Slow Quasi-Manifold

Three of Lighthill's many interests, (a) wave propagation in moving media, (b) acoustic streaming and related phenomena, and (c), in a very subtle and fascinating way, aerodynamic sound generation, are all turning out to be fundamental to understanding our atmospheric environment. Among other things there is the global-scale circulation that shapes the ozone layer and controls the rate of destruction of man-made chlorofluorocarbons (CFCs). Recent progress in this field is sketched, including progress in understanding the abstract structure of Hamiltonian theories of balanced motion, the so-called slow “manifold”. Here a generic phenomenon, “velocity splitting”, turns out to be intimately related to aerodynamic sound generation, particularly its generalization describing the spontaneous emission of inertia–gravity waves from unsteady vortical motions in stratified rotating flow. Received 5 January 1997 and accepted 2 May 1997     

Experimental study and direct numerical simulation of the evolution of disturbances in a viscous shock layer on a flat plate

The evolution of disturbances in a hypersonic viscous shock layer on a flat plate excited by slow-mode acoustic waves is considered numerically and experimentally. The parameters measured in the experiments performed with a free-stream Mach number M _∞ = 21 and Reynolds number Re _L = 1.44 · 10⁵ are the transverse profiles of the mean density and Mach number, the spectra of density fluctuations, and growth rates of natural disturbances. Direct numerical simulation of propagation of disturbances is performed by solving the Navier-Stokes equations with a high-order shock-capturing scheme. The numerical and experimental data characterizing the mean flow field, intensity of density fluctuations, and their growth rates are found to be in good agreement. Possible mechanisms of disturbance generation and evolution in the shock layer at hypersonic velocities are discussed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 3–15, September–October, 2006.     

Application of an acoustic analogy to PIV data from rectangular cavity flows

The present paper describes a method to derive information about the acoustic emission of a flow using particle image velocimetry (PIV) data. The advantage of the method is that it allows studying sound sources, the related flow phenomena and their acoustic radiation into the far field, simultaneously. In a first step the time history of two-dimensional instantaneous pressure fields is derived from planar PIV data. In a successive step the Curle’s acoustic analogy is applied to the pressure data to obtain the acoustic radiation of the flow. The test cases studied here are two rectangular cavity flows at very low Mach number with different aspect ratios L/H. The main sound source is located at the cavity trailing edge and it is due to the impingement of vortices shed in the shear layer. It is shown that the flow emits sound with a main directivity in the upstream direction for the smaller aspect ratio and the directivity is more uniform for the larger aspect ratio. In the latter case the acoustic pressure spectra has a broader character due to the impact of the downstream recirculation zone onto the shear layer instabilities, destroying their regular pattern and alternating the main sound source.     

A Three-Dimensional Hybrid LES-Acoustic Analogy Method for Predicting Open-Cavity Noise

A three-dimensional (3D) hybrid LES-acoustic analogy method for computational aeroacoustics (CAA) is presented for the prediction of open-cavity noise. The method uses large-eddy simulation (LES) to compute the acoustic source while the Ffowcs Williams-Hawkings (FW-H) acoustic analogy is employed for the prediction of the far-field sound. As a comparison, a two-dimensional (2D) FW-H analogy is also included. The hybrid method has been assessed in an open-cavity flow at a Mach number of 0.85 and a Reynolds number of Re=1.36×10⁶, where some experimental data are available for comparison. The study has identified some important technical issues in the application of the FW-H acoustic analogy to cavity noise prediction and CAA in general, including the proper selection of the integration period and the modes of sound sources in the frequency domain. The different nature of 2D and 3D wave propagation is also highlighted, which calls for a matching acoustic solver for each problem. The developed hybrid method has shown promise to be a feasible, accurate and computationally affordable approach for CAA.     

Numerical simulation of receptivity of a hypersonic boundary layer to acoustic disturbances

Direct numerical simulations of the evolution of disturbances in a viscous shock layer on a flat plate are performed for a free-stream Mach number M _∞ = 21 and Reynolds number Re _L = 1.44 · 10⁵. Unsteady Navier-Stokes equations are solved by a high-order shock-capturing scheme. Processes of receptivity and instability development in a shock layer excited by external acoustic waves are considered. Direct numerical simulations are demonstrated to agree well with results obtained by the locally parallel linear stability theory (with allowance for the shock-wave effect) and with experimental measurements in a hypersonic wind tunnel. Mechanisms of conversion of external disturbances to instability waves in a hypersonic shock layer are discussed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 84–91, May–June, 2007.     

Numerical simulation of flow‐induced noise using LES/SAS and Lighthill's acoustic analogy

We present a finite element (FE) formulation of Lighthill's acoustic analogy for the hybrid computation of noise generated by turbulent flows. In the present approach, the flow field is computed using large eddy simulation and scale adaptive simulation turbulence models. The acoustic propagation is obtained by solving the variational formulation of Lighthill's acoustic analogy with the FE method. In order to preserve the acoustic energy, we compute the inhomogeneous part of Lighthill's wave equation by applying the FE formulation on the fine flow grid. The resulting acoustic nodal loads are then conservatively interpolated to the coarser acoustic grid. Subsequently, the radiated acoustic field can be solved in both time and frequency domains. In the latter case, an enhanced perfectly matched layer technique is employed, allowing one to truncate the computational domain in the acoustic near field, without compromising the numerical solution. Our hybrid approach is validated by comparing the numerical results of the acoustic field induced by a corotating vortex pair with the corresponding analytical solution. To demonstrate the applicability of our scheme, we present full 3D numerical results for the computed acoustic field generated by the turbulent flow around square cylinder geometries. The sound pressure levels obtained compare well with measured values. Copyright © 2009 John Wiley & Sons, Ltd.     

Comparison of low Mach number models for natural convection problems

 We investigate in this paper two numerical methods for solving low Mach number compressible flows and their application to single-phase natural convection flow problems. The first method is based on an asymptotic model of the Navier–Stokes equations valid for small Mach numbers, whereas the second is based on the full compressible Navier–Stokes equations with particular care given to the discretization at low Mach numbers. These models are more general than the Boussinesq incompressible flow model, in the sense that they are valid even for cases in which the fluid is subjected to large temperature differences, that is when the compressibility of the fluid manifests itself through low Mach number effects. Numerical solutions are computed for a series of test problems with fixed Rayleigh number and increasing temperature differences, as well as for varying Rayleigh number for a given temperature difference. Numerical difficulties associated with low Mach number effects are discussed, as well as the accuracy of the approximations. Received on 17 January 2000     

Development of localized arc filament RF plasma actuators for high-speed and high Reynolds number flow control

Recently developed localized arc filament plasma actuators (LAFPAs) have shown tremendous control authority in high-speed and high Reynolds number flow for mixing enhancement and noise mitigation. Previously, these actuators were powered by a high-voltage pulsed DC plasma generator with low energy coupling efficiency of 5–10%. In the present work, a new custom-designed 8-channel pulsed radio frequency (RF) plasma generator has been developed to power up to 8 plasma actuators operated over a wide range of forcing frequencies (up to 50 kHz) and duty cycles (1–50%), and at high energy coupling efficiency (up to 80–85%). This reduces input electrical power requirements by approximately an order of magnitude, down to 12 W per actuator operating at 10% duty cycle. The new pulsed RF plasma generator is scalable to a system with a large number of channels. Performance of pulsed RF plasma actuators used for flow control was studied in a Mach 0.9 circular jet with a Reynolds number of about 623,000 and compared with that of pulsed DC actuators. Eight actuators were distributed uniformly on the perimeter of a 2.54-cm diameter circular nozzle extension. Both types of actuators coupled approximately the same amount of power to the flow, but with drastically different electrical inputs to the power supplies. Particle image velocimetry measurements showed that jet centerline Mach number decay produced by DC and RF actuators operating at the same forcing frequencies and duty cycles is very similar. At a forcing Strouhal number near 0.3, close to the jet column instability frequency, well-organized periodic structures, with similar patterns and dimensions, were generated in the jets forced by both DC and RF actuators. Far-field acoustic measurements demonstrated similar trends in the overall sound pressure level (OASPL) change produced by both types of actuators, resulting in OASPL reduction up to 1.2–1.5 dB in both cases. We conclude that pulsed RF actuators demonstrate flow control authority similar to pulsed DC actuators, with a significantly reduced power budget.     

Investigation of density ratio effects on normally injected cold jets into a hot cross flow

In this work, hydrodynamic effects of multiple squared cross-section cold jets inclined normally into a hot cross flow were computationally simulated using large eddy simulation (LES) approach. The finite volume method and the SIMPLE algorithm on a multi-block non-uniform staggered grid were applied. The jet into cross flow velocity ratio and the jet Reynolds number were 0.5 and 4,700, respectively. Simulations were performed for three different jets into cross flow temperature ratios, namely, 2, 1, and 0.5, corresponding to density ratios of 0.5, 1, and 2. The experimental and numerical results of Ajersch et al. (J. Turbomachinery 119(2), 330–342, 1997) at unit density ratio were used for validation purposes. The density ratio effects on the flow characteristics were investigated. It was shown that the density ratio has significant effects on the hydrodynamics of the jet into cross flow problems and even though the flow Mach number may be low, its effects due to temperature differences between the jet and the cross flow cannot be easily ignored.     

Large-eddy Simulation of Film Cooling Flows with Variable Density Jets

The present paper investigates the impact of the velocity and density ratio on the turbulent mixing process in gas turbine blade film cooling. A cooling fluid is injected from an inclined pipe at α=30° into a turbulent boundary layer profile at a freestream Reynolds number of Re_∞ = 400,000. This jet-in-a-crossflow (JICF) problem is investigated using large-eddy simulations (LES). The governing equations comprise the Navier–Stokes equations plus additional transport equations for several species to simulate a non-reacting gas mixture. A variation of the density ratio is simulated by the heat-mass transfer analogy, i.e., gases of different density are effused into an air crossflow at a constant temperature. An efficient large-eddy simulation method for low subsonic flows based on an implicit dual time-stepping scheme combined with low Mach number preconditioning is applied. The numerical results and experimental velocity data measured using two-component particle-image velocimetry (PIV) are in excellent agreement. The results show the dynamics of the flow field in the vicinity of the jet hole, i.e., the recirculation region and the inclination of the shear layers, to be mainly determined by the velocity ratio. However, evaluating the cooling efficiency downstream of the jet hole the mass flux ratio proves to be the dominant similarity parameter, i.e., the density ratio between the fluids and the velocity ratio have to be considered.     

A comparative study of sound generation by laminar,combusting and non-combusting jet flows

Sound production by two-dimensional, laminar jet flows with and without combustion is studied numerically and theoretically. The compressible Navier–Stokes, energy and progress variable equations are solved by resolving both the near field and the acoustics. The combusting jet flows are compared to non-combusting jets of the same jet Mach number, with the non-combusting, non-isothermal jets having the same steady temperature difference as the combusting jets. This infers that the magnitude of entropic and density disturbances is similar in some of the combusting and non-combusting cases. The flows are perturbed by a sinusoidal inlet velocity fluctuation at different Strouhal numbers. The computational domain is resolved to the far field in all cases, allowing direct examination of the sound radiated and its sources. Lighthill’s acoustic analogy is then solved numerically using Green’s functions. The radiated sound calculated using Lighthill’s equation is in good agreement with that from the simulations for all cases, validating the numerical solution of Lighthill’s equation. The contribution of the source terms in Dowling’s reformulation of Lighthill’s equation is then investigated. It is shown that the source term relating to changes in the momentum of density inhomogeneities is the dominant source term for all non-reacting, non-isothermal cases. Further, this source term has similar magnitude in the combusting cases and is one of the several source terms that have similar magnitude to the source term involving fluctuations in the heat release rate.     

平板湍流边界层的声辐射

基于平板湍流边界层的壁压起伏波数—频率谱 ,给出了一种湍流边界层声辐射的估算方法 ,并对光滑平板湍流边界层和平板表面粗糙度引起的湍流边界层声辐射进行了分析。结果表明 :湍流边界层声辐射是一种四极子声辐射 ,且其辐射声能集中于平板表面粗糙度引起的湍流边界层声辐射 ;光滑平板湍流边界层的声辐射也不可忽略。     

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].     

Application de méthodes intégrales au calcul du bruit de cavitéComputation of cavity noise using integral methods

A direct numerical simulation of the sound radiated by a flow over a 2-D subsonic cavity is performed using Computational AeroAcoustics tools. The simulation is consistent with Karamcheti's experiments. The numerical results are used as a reference to study two integral formulations: the Ffowcs Williams and Hawkings acoustic analogy and a wave extrapolation method based on FW-H equation. These hybrid approaches agree with the direct computation by DNS data and provide powerful tools to compute far-field noise. To cite this article: X. Gloerfelt et al., C. R. Mecanique 330 (2002) 13–20     

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.