Topology-based characterization of compressibility effects in mixing layers

2019 ◽  
Vol 874 ◽  
pp. 38-75 ◽  
Author(s):  
S. Arun ◽  
A. Sameen ◽  
B. Srinivasan ◽  
S. S. Girimaji

Direct numerical simulations of high-speed mixing layers are used to characterize the effects of compressibility on the basis of local streamline topology and vortical structure. Temporal simulations of the mixing layers are performed using a finite volume gas-kinetic scheme for convective Mach numbers ranging from $M_{c}=0.2$ to $M_{c}=1.2$. The focus of the study is on the transient development and the main objectives are to (i) investigate and characterize the turbulence suppression mechanism conditioned upon local streamline topology; and (ii) examine changes in the vortex vector field – distribution, magnitude and orientation – as a function of Mach number. We first reaffirm that kinetic energy suppression with increasing Mach number is due to a decrease in pressure–strain redistribution. Then, we examine the suppression mechanism conditioned upon topology and vortex structure. Conditional statistics indicate that (i) at a given Mach number, shear-dominated topologies generally exhibit more effective pressure–strain redistribution than vortical topologies; and (ii) for a given topology, the level of pressure–strain correlation mostly decreases with increasing Mach number. At each topology, with increasing Mach number, there is a corresponding decrease in turbulent shear stress and production leading to reduced kinetic energy. Further, as $M_{c}$ increases, the proportion of vortex-dominated regions in the flow increases, leading to further reduction in the turbulent kinetic energy of the flow. Then, the orientation of vortical structures and direction of fluid rotation are examined using the vortex vector approach of Tian et al. (J. Fluid Mech., vol. 849, 2018, pp. 312–339). At higher $M_{c}$, the vortex vectors tend to be more aligned in the streamwise direction in contrast to low $M_{c}$ wherein larger angles with streamwise direction are preferred. The connection between vortex orientation and kinetic energy production is also investigated. The findings lead to improved insight into turbulence suppression dynamics in high Mach number turbulent flows.

2014 ◽  
Vol 743 ◽  
pp. 554-584 ◽  
Author(s):  
J. O’Brien ◽  
J. Urzay ◽  
M. Ihme ◽  
P. Moin ◽  
A. Saghafian

AbstractThis study addresses the dynamics of backscatter of kinetic energy in the context of large-eddy simulations (LES) of high-speed turbulent reacting flows. A priori analyses of direct numerical simulations (DNS) of reacting and inert supersonic, time-developing, hydrogen–air turbulent mixing layers with complex chemistry and multicomponent diffusion are conducted here in order to examine the effects of compressibility and combustion on subgrid-scale (SGS) backscatter of kinetic energy. The main characteristics of the aerothermochemical field in the mixing layer are outlined. A selfsimilar period is identified in which some of the turbulent quantities grow in a quasi-linear manner. A differential filter is applied to the DNS flow field to extract filtered quantities of relevance for the large-scale kinetic-energy budget. Spatiotemporal analyses of the flow-field statistics in the selfsimilar regime are performed, which reveal the presence of considerable amounts of SGS backscatter. The dilatation field becomes spatially intermittent as a result of the high-speed compressibility effect. In addition, the large-scale pressure-dilatation work is observed to be an essential mechanism for the local conversion of thermal and kinetic energies. A joint probability density function (PDF) of SGS dissipation and large-scale pressure-dilatation work is provided, which shows that backscatter occurs primarily in regions undergoing volumetric expansion; this implies the existence of an underlying physical mechanism that enhances the reverse energy cascade. Furthermore, effects of SGS backscatter on the Boussinesq eddy viscosity are studied, and a regime diagram demonstrating the relationship between the different energy-conversion modes and the sign of the eddy viscosity is provided along with a detailed budget of the volume fraction in each mode. A joint PDF of SGS dissipation and SGS dynamic-pressure dilatation work is calculated, which shows that high-speed compressibility effects lead to a decorrelation between SGS backscatter and negative eddy viscosities, which increases for increasingly large values of the SGS Mach number and filter width. Finally, it is found that the combustion dynamics have a marginal impact on the backscatter and flow-dilatation distributions, which are mainly dominated by the high-Mach-number effects.


Fluids ◽  
2022 ◽  
Vol 7 (1) ◽  
pp. 34
Author(s):  
Hechmi Khlifi ◽  
Adnen Bourehla

This work focuses on the performance and validation of compressible turbulence models for the pressure-strain correlation. Considering the Launder Reece and Rodi (LRR) incompressible model for the pressure-strain correlation, Adumitroaie et al., Huang et al., and Marzougui et al., used different modeling approaches to develop turbulence models, taking into account compressibility effects for this term. Two numerical coefficients are dependent on the turbulent Mach number, and all of the remaining coefficients conserve the same values as in the original LRR model. The models do not correctly predict the compressible turbulence at a high-speed shear flow. So, the revision of these models is the major aim of this study. In the present work, the compressible model for the pressure-strain correlation developed by Khlifi−Lili, involving the turbulent Mach number, the gradient, and the convective Mach numbers, is used to modify the linear mean shear strain and the slow terms of the previous models. The models are tested in two compressible turbulent flows: homogeneous shear flow and the newly developed plane mixing layers. The predicted results of the proposed modifications of the Adumitroaie et al., Huang et al., and Marzougui et al., models and of its universal versions are compared with direct numerical simulation (DNS) and experiment data. The results show that the important parameters of compressibility in homogeneous shear flow and in the mixing layers are well predicted by the proposal models.


Author(s):  
Y. Jiang ◽  
N. Gurram ◽  
E. Romero ◽  
P. T. Ireland ◽  
L. di Mare

Slot film cooling is a popular choice for trailing edge cooling in high pressure (HP) turbine blades because it can provide more uniform film coverage compared to discrete film cooling holes. The slot geometry consists of a cut back in the blade pressure side connected through rectangular openings to the internal coolant feed passage. The numerical simulation of this kind of film cooling flows is challenging due to the presence of flow interactions like step flow separation, coolant-mainstream mixing and heat transfer. The geometry under consideration is a cutback surface at the trailing edge of a constant cross-section aerofoil. The cutback surface is divided into three sections separated by narrow lands. The experiments are conducted in a high speed cascade in Oxford Osney Thermo-Fluids Laboratory at Reynolds and Mach number distributions representative of engine conditions. The capability of CFD methods to capture these flow phenomena is investigated in this paper. The isentropic Mach number and film effectiveness are compared between CFD and pressure sensitive paint (PSP) data. Compared to steady k–ω SST method, Scale Adaptive Simulation (SAS) can agree better with the measurement. Furthermore, the profiles of kinetic energy, production and shear stress obtained by the steady and SAS methods are compared to identify the main source of inaccuracy in RANS simulations. The SAS method is better to capture the unsteady coolant-hot gas mixing and vortex shedding at the slot lip. The cross flow is found to affect the film significantly as it triggers flow separation near the lands and reduces the effectiveness. The film is non-symmetric with respect to the half-span plane and different flow features are present in each slot. The effect of mass flow ratio (MFR) on flow pattern and coolant distribution is also studied. The profiles of velocity, kinetic energy and production of turbulent energy are compared among the slots in detail. The MFR not only affects the magnitude but also changes the sign of production.


Author(s):  
D. Basu ◽  
A. Hamed ◽  
K. Das

This study deals with the computational grid requirements in multiscale simulations of separated turbulent flows at high Reynolds number. The two-equation k-ε based DES (Detached Eddy Simulation) model is implemented in a full 3-D Navier-Stokes solver and numerical results are presented for transonic flow solution over an open cavity. Results for the vorticity, pressure fluctuations, SPL (Sound Pressure level) spectra and for modeled and resolved TKE (Turbulent Kinetic Energy) are presented and compared with available experimental data and with LES results. The results indicate that grid resolution significantly influences the resolved scales and the peak amplitude of the unsteady sound pressure level (SPL) and turbulent kinetic energy spectra.


1997 ◽  
Vol 330 ◽  
pp. 307-338 ◽  
Author(s):  
A. SIMONE ◽  
G.N. COLEMAN ◽  
C. CAMBON

The influence of compressibility upon the structure of homogeneous sheared turbulence is investigated. For the case in which the rate of shear is much larger than the rate of nonlinear interactions of the turbulence, the modification caused by compressibility to the amplification of turbulent kinetic energy by the mean shear is found to be primarily reflected in pressure–strain correlations and related to the anisotropy of the Reynolds stress tensor, rather than in explicit dilatational terms such as the pressure–dilatation correlation or the dilatational dissipation. The central role of a ‘distortion Mach number’ Md =  S[lscr ]/a, where S is the mean strain or shear rate, [lscr ] a lengthscale of energetic structures, and a the sonic speed, is demonstrated. This parameter has appeared in previous rapid-distortion-theory (RDT) and direct-numerical-simulation (DNS) studies; in order to generalize the previous analyses, the quasi-isentropic compressible RDT equations are numerically solved for homogeneous turbulence subjected to spherical (isotropic) compression, one-dimensional (axial) compression and pure shear. For pure-shear flow at finite Mach number, the RDT results display qualitatively different behaviour at large and small non-dimensional times St: when St < 4 the kinetic energy growth rate increases as the distortion Mach number increases; for St > 4 the inverse occurs, which is consistent with the frequently observed tendency for compressibility to stabilize a turbulent shear flow. This ‘crossover’ behaviour, which is not present when the mean distortion is irrotational, is due to the kinematic distortion and the mean-shear-induced linear coupling of the dilatational and solenoidal fields. The relevance of the RDT is illustrated by comparison to the recent DNS results of Sarkar (1995), as well as new DNS data, both of which were obtained by solving the fully nonlinear compressible Navier–Stokes equations. The linear quasi-isentropic RDT and nonlinear non-isentropic DNS solutions are in good general agreement over a wide range of parameters; this agreement gives new insight into the stabilizing and destabilizing effects of compressibility, and reveals the extent to which linear processes are responsible for modifying the structure of compressible turbulence.


2008 ◽  
Vol 112 (1133) ◽  
pp. 381-394 ◽  
Author(s):  
G. M. Lilley

Abstract The present paper reviews and discusses the physical mechanisms of noise generation and reduction in turbulent flows with their applications towards aircraft noise reduction at takeoff and on the approach. This work began in 1948 when Lilley undertook an experimental investigation into the source of jet noise as a necessary precursor to finding methods for the reduction of high speed jet engine noise on civil jet airliners. Westley and Lilley completed this experimental programme in 1951, which included the design of a range of devices for high speed jet noise reduction. It was about this time that similar studies on jet noise were being started elsewhere and in particular by Lassiter and Hubbard in USA. The major contribution to the subject of turbulence as a source of noise came from Sir James Lighthill’s remarkable theory in 1952. In spite of the difficulties attached to theoretical and experimental studies on noise from turbulence, it is shown that with the accumulated knowledge on aerodynamic noise over the past 50 years, together with an optimisation of aircraft operations including flight trajectories, we are today on the threshold of approaching the design of commercial aircraft with turbofan propulsion engines that will not be heard above the background noise of the airport at takeoff and landing beyond 1-2km, from the airport boundary fence. It is evident that in the application of this work, which centres on the physical mechanisms relating to the generation of noise from turbulence and turbulent shear flows, to jet noise, there is not one unique mechanism of jet noise generation for all jet Mach numbers. This author in this publication has concentrated on what appears to be the dominant mechanism of noise generation from turbulence, where the mean convection speeds of the turbulence are subsonic. The noise generated at transonic and supersonic jet speeds invariably involves extra mechanisms, which are only briefly referred to here.


Author(s):  
T. Kawamura ◽  
T. Nakatani

Direct numerical simulations of homogeneous shear turbulent flows containing deformable bubbles were carried out for clarifying the mechanism of drag reduction by microbubbles. The results show that presence of bubbles can suppress or enhance the development of turbulence depending on condition. The dissipation rate of turbulent kinetic energy is always increased by bubbles, while the production rate can be either increased or decreased depending on the turbulent and shear Reynolds numbers. As a result, the growth rate of turbulent kinetic energy can be either increased or decreased by bubbles depending on conditions. It was shown that the production rate tends to decrease at smaller shear Reynolds number, at larger turbulent Reynolds number, and at larger Weber number. Based on the results, a hypothetical model to explain the dependency on the Reynolds numbers has been proposed.


2019 ◽  
Vol 91 (7) ◽  
pp. 967-976
Author(s):  
Chang-Fei Zhuo ◽  
Ming-Xiao Wang ◽  
Wen-Jin Yao ◽  
Wen-ke Xu

Purpose The purpose of this paper is to study the operation performance of the high-speed ramjet kinetic energy projectile using solid fuel ramjet as power plant that is a new short-range and small caliber projectile. Design/methodology/approach The numerical investigation on combustion characteristic of polyethylene in high-speed ramjet kinetic energy projectile is carried out in this paper. The flow characteristics’ differences are analyzed when ramjet works or do not work, and both the combustion characteristics and propulsive performance are analyzed when ramjet works. Findings The results show that with the increase of the abscissa x, the flame front is close to solid fuel surface at first and then keeps away from solid fuel surface. With the increase of the abscissa x, the temperature of solid fuel surface and regression rate of solid fuel continues to increase before re-attachment point and then decreases, which a maximum locate at the re-attachment point. Both the average temperature and the regression rate on the surface of the solid fuel tend to rise as the increase of inflow Mach number. As the inflow of Mach number increases, the mass flow rate of gaseous fuel increases. Practical implications The research results can provide useful database for the subsequent research on high-speed ramjet kinetic energy projectile. Originality/value This paper studies the operation characteristics of the ramjet projectile, especially the effect of the change of the flight velocity on the performance of high-speed ramjet projectile.


Author(s):  
Takeo Kajishima ◽  
Takashi Ohta

Flow field of low Mach number (e.g. M &lt;0.3) is usually simulated by the incompressible flow scheme due to the severe limitation of time-increment in the compressible flow scheme. In this work, we propose a modification to the usual incompressible scheme, based on the elliptic equation for pressure, to improve the accuracy for turbulent flows considering weak compressibility. Two examples will be shown to validate our method. (1) LES (Large-Eddy Simulation) was conducted for turbulent flow around NACA0012 airfoil. Particular attention was focused on the influence of compressibility, despite the low Mach number range. In addition, new subgrid scale model of one-equation type using dynamic procedure was compared with traditional Smagorinsky model. Our method successfully reproduced the separation bubble near the leading edge, resulted in the improvement in the intensity of pressure fluctuation. (2) DNS (Direct Numerical Simulation) of turbulent flow in a plane channel is carried out, taking wall temperature difference into account. As a result of the density fluctuation in near-wall eddies, asymmetric profiles are observed in turbulence statistics. By the 4-quadrant analysis of turbulent shear stress, it is found that the ejection events in the vicinity of the walls are particularly affected by the density variation.


1972 ◽  
Vol 54 (1) ◽  
pp. 39-48 ◽  
Author(s):  
James M. Wallace ◽  
Helmut Eckelmann ◽  
Robert S. Brodkey

Hot-film measurements in a fully developed channel flow have been made in an attempt to gain more insight into the process of Reynolds stress production. The background for this effort is the observation of a certain sequence of events (deceleration, ejection and sweep) in the wall region of turbulent flows by Corino (1965) and Corino & Brodkey (1969). The instantaneous product signal uv was classified according to the sign of its components u and v, and these classified portions were then averaged to obtain their contributions to the Reynolds stress $-\rho\overline{uv} $. The signal was classified into four categories; the two main ones were that with u negative and v positive, which can be associated with the ejection-type motion of Corino & Brodkey (1969), and that with u positive and v negative, associated with the sweep-type motion. It was found that over the wall region investigated, 3·5 [les ] y [les ] 100, these two types of motion give rise to a stress considerably greater than the total Reynolds stress. Two other types of motion, (i) u negative, v negative, corresponding to low-speed fluid deflected towards the wall, and (ii) u positive, v positive, corresponding to high-speed fluid reflected outwards from the wall, were found to account for the ‘excess’ stress produced by the first two categories, which give contributions of opposite sign.The autocorrelations of the classified portions of uv were obtained to determine the relative time scales of these four types of motion. The positive stress producing motions (u < 0, v > 0 and u > 0, v < 0) were found to have significantly larger time scales than the negative stress producing motions (u < 0, v < 0 and u > 0, v > 0). It was further surmised that turbulent energy dissipation is associated with the Reynolds stress producing motions, since these result in localized shear regions in which the dissipation is several orders of magnitude greater than the average dissipation at the wall.


Sign in / Sign up

Export Citation Format

Share Document