Direct numerical simulation of compressible turbulence in a homogeneous shear flow

Author(s):  
S. Sarkar ◽  
G. Erlebacher
1993 ◽  
Vol 256 ◽  
pp. 443-485 ◽  
Author(s):  
G. A. Blaisdell ◽  
N. N. Mansour ◽  
W. C. Reynolds

Compressibility effects within decaying isotropic turbulence and homogeneous turbulent shear flow have been studied using direct numerical simulation. The objective of this work is to increase our understanding of compressible turbulence and to aid the development of turbulence models for compressible flows. The numerical simulations of compressible isotropic turbulence show that compressibility effects are highly dependent on the initial conditions. The shear flow simulations, on the other hand, show that measures of compressibility evolve to become independent of their initial values and are parameterized by the root mean square Mach number. The growth rate of the turbulence in compressible homogeneous shear flow is reduced compared to that in the incompressible case. The reduced growth rate is the result of an increase in the dissipation rate and energy transfer to internal energy by the pressure–dilatation correlation. Examination of the structure of compressible homogeneous shear flow reveals the presence of eddy shocklets, which are important for the increased dissipation rate of compressible turbulence.


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.


2017 ◽  
Vol 827 ◽  
pp. 225-249 ◽  
Author(s):  
Atsushi Sekimoto ◽  
Javier Jiménez

Unstable equilibrium solutions in a homogeneous shear flow with sinuous (streamwise-shift-reflection and spanwise-shift-rotation) symmetry are numerically found in large-eddy simulations (LES) with no kinetic viscosity. The small-scale properties are determined by the mixing length scale $l_{S}$ used to define eddy viscosity, and the large-scale motion is induced by the mean shear at the integral scale, which is limited by the spanwise box dimension $L_{z}$. The fraction $R_{S}=L_{z}/l_{S}$, which plays the role of a Reynolds number, is used as a numerical continuation parameter. It is shown that equilibrium solutions appear by a saddle-node bifurcation as $R_{S}$ increases, and that the flow structures resemble those in plane Couette flow with the same sinuous symmetry. The vortical structures of both lower- and upper-branch solutions become spontaneously localised in the vertical direction. The lower-branch solution is an edge state at low $R_{S}$, and takes the form of a thin critical layer as $R_{S}$ increases, as in the asymptotic theory of generic shear flow at high Reynolds numbers. On the other hand, the upper-branch solutions are characterised by a tall velocity streak with multiscale multiple vortical structures. At the higher end of $R_{S}$, an incipient multiscale structure is found. The LES turbulence occasionally visits vertically localised states whose vortical structure resembles the present vertically localised LES equilibria.


2009 ◽  
Vol 629 ◽  
pp. 25-39 ◽  
Author(s):  
P. GUALTIERI ◽  
F. PICANO ◽  
C. M. CASCIOLA

Recently, clustering of inertial particles in turbulence has been thoroughly analysed for statistically homogeneous isotropic flows. Phenomenologically, spatial homogeneity of particle configurations is broken by the advection of a range of eddies determined by the Stokes relaxation time of the particles. This in turn results in a multi-scale distribution of local particle concentration and voids. Much less is known concerning anisotropic flows. Here, by addressing direct numerical simulations (DNS) of a statistically steady particle-laden homogeneous shear flow, we provide evidence that the mean shear preferentially orients particle patterns. By imprinting anisotropy on large-scale velocity fluctuations, the shear indirectly affects the geometry of the clusters. Quantitative evaluation is provided by a purposely designed tool, the angular distribution function (ADF) of particle pairs, which allows to address the anisotropy content of particle aggregates on a scale-by-scale basis. The data provide evidence that, depending on the Stokes relaxation time of the particles, anisotropic clustering may occur even in the range of scales in which the carrier phase velocity field is already recovering isotropy. The strength of the singularity in the anisotropic component of the ADF quantifies the level of fine-scale anisotropy, which may even reach values of more than 30% direction-dependent variation in the probability to find two closeby particles at viscous-scale separation.


2001 ◽  
Author(s):  
Shriram B. Pillapakkam ◽  
Pushpendra Singh

Abstract A three dimensional finite element scheme for Direct Numerical Simulation (DNS) of viscoelastic two phase flows is implemented. The scheme uses the Level Set Method to track the interface and the Marchuk-Yanenko operator splitting technique to decouple the difficulties associated with the governing equations. Using this numerical scheme, the shape of Newtonian drops in a simple shear flow of viscoelastic fluid and vice versa are analyzed as a function of Capillary number, Deborah number and polymer concentration. The viscoelastic fluid is modeled via the Oldroyd-B model. The role of viscoelastic stresses in deformation of a drop subjected to simple shear flow and its effect on the steady state shape is analyzed. Our results compare favorably with existing experimental data and also help in understanding the role of viscoelastic stresses in drop deformation.


Sign in / Sign up

Export Citation Format

Share Document