scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Numerical simulation of compressible homogeneous flows in the turbulent regime

01 Sep 1987-Journal of Fluid Mechanics (Cambridge University Press)-Vol. 181, Iss: -1, pp 441-466
TL;DR: In this paper, a pseudospectral simulation of turbulent homogeneous flows with r.m. velocities of the order of the speed of sound was performed using the Navier-Stokes equations.
Abstract: Compressible flows with r.m.8. velocities of the order of the speed of sound are studied with direct numerical simulations using a pseudospectral method. We concentrate on turbulent homogeneous flows in the two-dimensional case. The fluid obeys the Navier-Stokes equations for a perfect gas, and viscous terms are included explicitly. No modelling of small scales is used. We show that the behaviour of the flow differs sharply at low compared with high r.m.9. Mach number Ma, with a transition at Mu = 0.3. In the large scales, temporal exchanges between longitudinal and solenoidal modes of energy retain an acoustical character; they lead to a slowing down of the decrease of the Mach number with time, which occurs with interspersed plateaux corresponding to quiescent periods. When the flow is initially supersonic, the small scales are dominated by shocks behind which vortices form. This vortex production is particularly prominent, when two strong shocks collide, with the onset of shear turbulence in the region downstream of the collision. However, at the resolutions reached by our code on a 256 x 256 uniform grid, this mechanism proves insufficient to bring vortices into equipartition with shocks in the small-scale tail of the energy spectrum.
Citations
More filters
Journal ArticleDOI
TL;DR: In this paper, a compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor.
Abstract: New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

714 citations

Journal ArticleDOI
TL;DR: In this article, the authors derive a shock capturing tool able to treat turbulence with minimum dissipation out of the shock for a large-eddy simulation (LES) of the interaction.

605 citations

Journal ArticleDOI
TL;DR: In this article, the authors report the results of numerical experiments in two opposite regimes: A ~ 1 and A 1, where A is the initial Alfvenic Mach number, the ratio of the rms velocity to the Alfven speed.
Abstract: Supersonic random motions are observed in dark clouds and are traditionally interpreted as Alfven waves, but the possibility that these motions are super-Alfvenic has not been ruled out. In this work we report the results of numerical experiments in two opposite regimes: A ~ 1 and A 1, where A is the initial Alfvenic Mach number—the ratio of the rms velocity to the Alfven speed. Our results show that models with A 1 are consistent with the observed properties of molecular clouds that we have tested (statistics of extinction measurements, distribution of integrated antenna temperature, Zeeman-splitting measurements of magnetic field strength, line width versus integrated antenna temperature of molecular emission-line spectra, statistical B-n relation, and scatter in that relation), while models withA ~ 1 have properties that are in conflict with the observations. We find that both the density and the magnetic field in molecular clouds may be very intermittent. The statistical distributions of the magnetic field and gas density are related by a power law, with an index that decreases with time in experiments with decaying turbulence. After about one dynamical time it stabilizes at B ∝ n0.4. Magnetically dominated cores form early in the evolution, while later on the intermittency in the density field wins out, and also cores with a weak field can be generated by mass accretion along magnetic field lines.

425 citations


Cites methods from "Numerical simulation of compressibl..."

  • ...D) simulations of turbulence with rms Mach number larger than one were performed by Passot & Pouquet (1987)....

    [...]

Journal ArticleDOI
TL;DR: In this paper, a dilatation dissipation model for high Reynolds number compressible turbulence is introduced, which is predicated on the existence of shocklike structures embedded within energetic turbulent eddies.
Abstract: In this paper a concept of dilatation dissipation ed for high Reynolds number compressible turbulence is introduced. The concept is predicated on the existence of shocklike structures embedded within energetic turbulent eddies. A parametric expression for ed is found that contains calculable parameters of a turbulent field: turbulence energy and length scale, rms (turbulent) Mach number, and the kurtosis of the fluctuating velocity. The dilatation dissipation is incorporated in a second‐order closure model for compressible mixing layers and model predictions of mean and turbulence quantities are presented and, where possible, compared with experiments. It is shown that the model is capable of predicting the reduction of layer growth rates as a function of the convective Mach number Mc in accordance with Papamoschou–Roshko experiments; the computations are also shown to compare well with available measurements of Reynolds stresses at Mc=0.5–0.86. Finally, the physical implications of the new model and resu...

407 citations

References
More filters
Journal ArticleDOI
TL;DR: In this paper, a theory for estimating the sound radiated from a fluid flow, with rigid boundaries, which as a result of instability contains regular fluctuations or turbulence is initiated, based on the equations of motion of a gas.
Abstract: A theory is initiated, based on the equations of motion of a gas, for the purpose of estimating the sound radiated from a fluid flow, with rigid boundaries, which as a result of instability contains regular fluctuations or turbulence. The sound field is that which would be produced by a static distribution of acoustic quadrupoles whose instantaneous strength per unit volume is ρv i v j + p ij - a 2 0 ρ δ ij , where ρ is the density, v i the velocity vector, p ij the compressive stress tensor, and a 0 the velocity of sound outside the flow. This quadrupole strength density may be approximated in many cases as ρ 0 v i v j . The radiation field is deduced by means of retarded potential solutions. In it, the intensity depends crucially on the frequency as well as on the strength of the quadrupoles, and as a result increases in proportion to a high power, near the eighth, of a typical velocity U in the flow. Physically, the mechanism of conversion of energy from kinetic to acoustic is based on fluctuations in the flow of momentum across fixed surfaces, and it is explained in § 2 how this accounts both for the relative inefficiency of the process and for the increase of efficiency with U . It is shown in § 7 how the efficiency is also increased, particularly for the sound emitted forwards, in the case of fluctuations convected at a not negligible Mach number.

4,697 citations

Journal ArticleDOI
TL;DR: In this paper, it was shown that two-dimensional turbulence has both kinetic energy and mean square vorticity as inviscid constants of motion, and two formal inertial ranges, E(k)∼e2/3k−5/3/3, where e is the rate of cascade of kinetic energy per unit mass, η is the time taken to reach a cascade of mean square velocity, and k is the kinetic energy of the entire mass.
Abstract: Two‐dimensional turbulence has both kinetic energy and mean‐square vorticity as inviscid constants of motion. Consequently it admits two formal inertial ranges, E(k)∼e2/3k−5/3 and E(k)∼η2/3k−3, where e is the rate of cascade of kinetic energy per unit mass, η is the rate of cascade of mean‐square vorticity, and the kinetic energy per unit mass is ∫0∞E(k) dk. The −53 range is found to entail backward energy cascade, from higher to lower wavenumbers k, together with zero‐vorticity flow. The −3 range gives an upward vorticity flow and zero‐energy flow. The paradox in these results is resolved by the irreducibly triangular nature of the elementary wavenumber interactions. The formal −3 range gives a nonlocal cascade and consequently must be modified by logarithmic factors. If energy is fed in at a constant rate to a band of wavenumbers ∼ki and the Reynolds number is large, it is conjectured that a quasi‐steady‐state results with a −53 range for k « ki and a −3 range for k » ki, up to the viscous cutoff. The t...

2,950 citations

Journal ArticleDOI
TL;DR: The theory of sound generated aerodynamically is extended by taking into account the statistical properties of turbulent airflows, from which the sound radiated (without the help of solid boundaries) is called aerodynamic noise as mentioned in this paper.
Abstract: The theory of sound generated aerodynamically is extended by taking into account the statistical properties of turbulent airflows, from which the sound radiated (without the help of solid boundaries) is called aerodynamic noise. The theory is developed with special reference to the noise of jets, for which a detailed comparison with experiment is made (§7 for subsonic jets, §8 for supersonic ones). The quadrupole distribution of part I (Lighthill 1952) is shown to behave (see §3) as if it were concentrated into independent point quadrupoles, one in each ‘average eddy volume’. The sound field of each of these is distorted, in favour of downstream emission, by the general downstream motion of the eddy, in accordance with the quadrupole convection theory of part I. This explains, for jet noise, the marked preference for downstream emission, and its increase with jet velocity. For jet velocities considerably greater than the atmospheric speed of sound, the ‘Mach number of convection’ M c may exceed I in parts of the jet, and then the directional maximum for emission from these parts of the jet is at an angle of sec -1 ( M c ) to the axis (§8). Although turbulence without any mean flow has an acoustic power output, which was calculated to a rough approximation from the expressions of part I by Proudman (1952) (see also § 4 below), nevertheless, turbulence of given intensity can generate more sound in the presence of a large mean shear (§ 5). This sound has a directional maximum at 45° (or slightly less, due to the quadrupole convection effect) to the shear layer. These results follow from the fact that the most important term in the rate of change of momentum flux is the product of the pressure and the rate of strain (see figure 2). The higher frequency sound from the heavily sheared mixing region close to the orifice of a jet is found to be of this character. But the lower frequency sound from the fully turbulent core of the jet, farther downstream, can be estimated satisfactorily (§7) from Proudman’s results, which are here reinterpreted (§5) in terms of sound generated from combined fluctuations of pressure and rate of shear in the turbulence. The acoustic efficiency of the jet is of the order of magnitude 10 -4 M 5 , where M is the orifice Mach number. However, the good agreement, as regards total acoustic power output, with the dimensional considerations of part I, is partly fortuitous. The quadrupole convection effect should produce an increase in the dependence of acoustic power on the jet velocity above the predicted U 8 law. The experiments show that (largely cancelling this) some other dependence on velocity is present, tending to reduce the intensity, at the stations where the convection effect would be absent, below the U 8 law. At these stations (at 90° to the jet) proportionality to about U 6.5 is more common. A suggested explanation of this, compatible with the existing evidence, is that at higher Mach numbers there may be less turbulence (especially for larger values of nd / U , where n is frequency and d diameter), because in the mixing region, where the turbulence builds up, it is losing energy by sound radiation. This would explain also the slow rate of spread of supersonic mixing regions, and, indeed, is not incompatible with existing rough explanations of that phenomenon. A consideration (§6) of whether the terms other than momentum flux in the quadrupole strength density might become important in heated jets indicates that they should hardly ever be dominant. Accordingly, the physical explanation (part I) of aerodynamic sound generation still stands. It is re-emphasized, however, that whenever there is a fluctuating force between the fluid and a solid boundary, a dipole radiation will result which may be more efficient than the quadrupole radiation, at least at low Mach numbers.

1,479 citations