Evolution of enstrophy in shock/homogeneous turbulence interaction

2012 ◽  
Vol 707 ◽  
pp. 74-110 ◽  
Author(s):  
Krishnendu Sinha
Keyword(s):  
Shock Wave ◽  
Numerical Simulations ◽  
High Speed ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Interaction Analysis ◽  
Direct Numerical Simulations ◽  
Homogeneous Isotropic Turbulence ◽  
Mean Flow ◽  
Linear Interaction

AbstractInteraction of turbulent fluctuations with a shock wave plays an important role in many high-speed flow applications. This paper studies the amplification of enstrophy, defined as mean-square fluctuating vorticity, in homogeneous isotropic turbulence passing through a normal shock. Linearized Navier–Stokes equations written in a frame of reference attached to the unsteady shock wave are used to derive transport equations for the vorticity components. These are combined to obtain an equation that describes the evolution of enstrophy across a time-averaged shock wave. A budget of the enstrophy equation computed using results from linear interaction analysis and data from direct numerical simulations identifies the dominant physical mechanisms in the flow. Production due to mean flow compression and baroclinic torques are found to be the major contributors to the enstrophy amplification. Closure approximations are proposed for the unclosed correlations in the production and baroclinic source terms. The resulting model equation is integrated to obtain the enstrophy jump across a shock for a range of upstream Mach numbers. The model predictions are compared with linear theory results for varying levels of vortical and entropic fluctuations in the upstream flow. The enstrophy model is then cast in the form of$k$–$\epsilon $equations and used to compute the interaction of homogeneous isotropic turbulence with normal shocks. The results are compared with available data from direct numerical simulations. The equations are further used to propose a model for the amplification of turbulent viscosity across a shock, which is then applied to a canonical shock–boundary layer interaction. It is shown that the current model is a significant improvement over existing models, both for homogeneous isotropic turbulence and in the case of complex high-speed flows with shock waves.

Turbulent energy flux generated by shock/homogeneous-turbulence interaction

2016 ◽  
Vol 796 ◽  
pp. 113-157 ◽  
Author(s):  
Russell Quadros ◽  
Krishnendu Sinha ◽  
Johan Larsson
Keyword(s):  
Shock Wave ◽  
Turbulent Flows ◽  
Energy Flux ◽  
High Speed ◽  
Isotropic Turbulence ◽  
Interaction Analysis ◽  
Turbulent Energy ◽  
Angle Of Incidence ◽  
Two Dimensional ◽  
Linear Interaction

High-speed turbulent flows with shock waves are characterized by high localized surface heat transfer rates. Computational predictions are often inaccurate due to the limitations in modelling of the unclosed turbulent energy flux in the highly non-equilibrium regions of shock interaction. In this paper, we investigate the turbulent energy flux generated when homogeneous isotropic turbulence passes through a nominally normal shock wave. We use linear interaction analysis where the incoming turbulence is idealized as being composed of a collection of two-dimensional planar vorticity waves, and the shock wave is taken to be a discontinuity. The nature of the postshock turbulent energy flux is predicted to be strongly dependent on the angle of incidence of the incoming waves. The energy flux correlation is also decomposed into its vortical, entropy and acoustic contributions to understand its rapid non-monotonic variation behind the shock. Three-dimensional statistics, calculated by integrating two-dimensional results over a prescribed upstream energy spectrum, are compared with available data from direct numerical simulations. A detailed budget of the governing equation is also considered in order to gain insight into the underlying physics.

Turbulence structure behind the shock in canonical shock–vortical turbulence interaction

2014 ◽  
Vol 756 ◽  
Author(s):  
Jaiyoung Ryu ◽  
Daniel Livescu
Keyword(s):  
Shock Wave ◽  
Mach Number ◽  
Isotropic Turbulence ◽  
Interaction Analysis ◽  
Longitudinal Velocity ◽  
Wave Structure ◽  
Turbulence Structure ◽  
Turbulent Structures ◽  
Linear Interaction ◽  
Third Invariant

AbstractThe interaction between vortical isotropic turbulence (IT) and a normal shock wave is studied using direct numerical simulation (DNS) and linear interaction analysis (LIA). In previous studies, agreement between the simulation results and the LIA predictions has been limited and, thus, the significance of LIA has been underestimated. In this paper, we present high-resolution simulations which accurately solve all flow scales (including the shock-wave structure) and extensively cover the parameter space (the shock Mach number, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}M_s$, ranges from 1.1 to 2.2 and the Taylor Reynolds number, ${\mathit{Re}}_{\lambda }$, ranges from 10 to 45). The results show, for the first time, that the turbulence quantities from DNS converge to the LIA solutions as the turbulent Mach number, $M_t$, becomes small, even at low upstream Reynolds numbers. The classical LIA formulae are extended to compute the complete post-shock flow fields using an IT database. The solutions, consistent with the DNS results, show that the shock wave significantly changes the topology of the turbulent structures, with a symmetrization of the third invariant of the velocity gradient tensor and ($M_s$-mediated) of the probability density function (PDF) of the longitudinal velocity derivatives, and an $M_s$-dependent increase in the correlation between strain and rotation.

The analysis and modelling of dilatational terms in compressible turbulence

1991 ◽  
Vol 227 ◽  
pp. 473-493 ◽  
Author(s):  
S. Sarkar ◽  
G. Erlebacher ◽  
M. Y. Hussaini ◽  
H. O. Kreiss
Keyword(s):  
Asymptotic Analysis ◽  
Numerical Simulations ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Mixing Layer ◽  
Direct Numerical Simulations ◽  
Homogeneous Turbulence ◽  
Compressible Turbulence ◽  
Quasi Equilibrium ◽  
Compressible Mixing Layer

It is shown that the dilatational terms that need to be modelled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of the compressible velocity field, which generates these dilatational terms, is explored by asymptotic analysis of the compressible Navier-Stokes equations in the case of homogeneous turbulence. The lowest-order equations for the compressible field are solved and explicit expressions for some of the associated one-point moments are obtained. For low Mach numbers, the compressible mode has a fast timescale relative to the incompressible mode. Therefore, it is proposed that, in moderate Mach number homogeneous turbulence, the compressible component of the turbulence is in quasi-equilibrium with respect to the incompressible turbulence. A non-dimensional parameter which characterizes this equilibrium structure of the compressible mode is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.

Direct Numerical Simulations of Colliding Particles Suspended in Homogeneous Isotropic Turbulence

2009 ◽  
Author(s):  
M. Ernst ◽  
M. Sommerfeld
Keyword(s):  
Numerical Simulations ◽  
Isotropic Turbulence ◽  
Three Dimensional ◽  
Stokes Number ◽  
Direct Numerical Simulations ◽  
Spherical Particles ◽  
Collision Rate ◽  
Homogeneous Isotropic Turbulence ◽  
And Performance ◽  
Fluid Behaviour

Direct numerical simulations of particle-laden homogeneous isotropic turbulence are performed to characterize the collision rate as a function of different particle properties. The fluid behaviour is computed using a three-dimensional Lattice Boltzmann Method including a spectral forcing scheme to generate the turbulence field. Under assumption of mass points, the transport of spherical particles is modelled in a Lagrangian frame of reference. In the simulations the influence of the particle phase on the fluid flow is neglected. The detection and performance of inelastic interparticle collisions are based on a deterministic collision model. Different studies with monodisperse particles are considered. According to the executed simulations, particles with small Stokes number possess a collision rate similar to the prediction of Saffman and Turner [1], whereas particles with larger Stokes numbers behave similarly to the theory of Abrahamson [2].

Sign in / Sign up

Export Citation Format

Share Document