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.

Settling velocity and concentration distribution of heavy particles in homogeneous isotropic turbulence

1993 ◽  
Vol 256 ◽  
pp. 27-68 ◽  
Author(s):  
Lian-Ping Wang ◽  
Martin R. Maxey
Keyword(s):  
Numerical Simulations ◽  
Drag Force ◽  
Settling Velocity ◽  
Isotropic Turbulence ◽  
Concentration Field ◽  
Terminal Velocity ◽  
Direct Numerical Simulations ◽  
Flow Dynamics ◽  
Homogeneous Turbulence ◽  
Heavy Particles

The average settling velocity in homogeneous turbulence of a small rigid spherical particle, subject to a Stokes drag force, has been shown to differ from that in still fluid owing to a bias from the particle inertia (Maxey 1987). Previous numerical results for particles in a random flow field, where the flow dynamics were not considered, showed an increase in the average settling velocity. Direct numerical simulations of the motion of heavy particles in isotropic homogeneous turbulence have been performed where the flow dynamics are included. These show that a significant increase in the average settling velocity can occur for particles with inertial response time and still-fluid terminal velocity comparable to the Kolmogorov scales of the turbulence. This increase may be as much as 50% of the terminal velocity, which is much larger than was previously found. The concentration field of the heavy particles, obtained from direct numerical simulations, shows the importance of the inertial bias with particles tending to collect in elongated sheets on the peripheries of local vortical structures. This is coupled then to a preferential sweeping of the particles in downward moving fluid. Again the importance of Kolmogorov scaling to these processes is demonstrated. Finally, some consideration is given to larger particles that are subject to a nonlinear drag force where it is found that the nonlinearity reduces the net increase in settling velocity.

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.

scholarly journals Triple Decomposition of Velocity Gradient Tensor in Compressible Turbulence

Fluids ◽  
2021 ◽  
Vol 6 (3) ◽  
pp. 98
Author(s):  
Radouan Boukharfane ◽  
Aimad Er-raiy ◽  
Linda Alzaben ◽  
Matteo Parsani
Keyword(s):  
Rigid Body ◽  
Numerical Simulations ◽  
Scaling Laws ◽  
Isotropic Turbulence ◽  
Three Dimensional ◽  
Direct Numerical Simulations ◽  
Compressible Turbulence ◽  
Local Motion ◽  
Eighth Order ◽  
Compressibility Effects

The decomposition of the local motion of a fluid into straining, shearing, and rigid-body rotation is examined in this work for a compressible isotropic turbulence by means of direct numerical simulations. The triple decomposition is closely associated with a basic reference frame (BRF), in which the extraction of the biasing effect of shear is maximized. In this study, a new computational and inexpensive procedure is proposed to identify the BRF for a three-dimensional flow field. In addition, the influence of compressibility effects on some statistical properties of the turbulent structures is addressed. The direct numerical simulations are carried out with a Reynolds number that is based on the Taylor micro-scale of Reλ=100 for various turbulent Mach numbers that range from Mat=0.12 to Mat=0.89. The DNS database is generated with an improved seventh-order accurate weighted essentially non-oscillatory scheme to discretize the non-linear advective terms, and an eighth-order accurate centered finite difference scheme is retained for the diffusive terms. One of the major findings of this analysis is that regions featuring strong rigid-body rotations or straining motions are highly spatially intermittent, while most of the flow regions exhibit moderately strong shearing motions in the absence of rigid-body rotations and straining motions. The majority of compressibility effects can be estimated if the scaling laws in the case of compressible turbulence are rescaled by only considering the solenoidal contributions.

scholarly journals Direct numerical simulations of laminar and transitional flows in diverging pipes

2019 ◽  
Vol 30 (1) ◽  
pp. 75-92 ◽  
Author(s):  
Dhanush Vittal Shenoy ◽  
Mostafa Safdari Shadloo ◽  
Jorge Peixinho ◽  
Abdellah Hadjadj
Keyword(s):  
Velocity Profile ◽  
Numerical Simulations ◽  
Stokes Equations ◽  
Spectral Element ◽  
Skin Friction Coefficient ◽  
Finite Amplitude ◽  
Direct Numerical Simulations ◽  
Recirculation Region ◽  
Cross Sectional ◽  
Content Type

Purpose Fluid flows in pipes whose cross-sectional area are increasing in the stream-wise direction are prone to separation of the recirculation region. This paper aims to investigate such fluid flow in expansion pipe systems using direct numerical simulations. The flow in circular diverging pipes with different diverging half angles, namely, 45, 26, 14, 7.2 and 4.7 degrees, are considered. The flow is fed by a fully developed laminar parabolic velocity profile at its inlet and is connected to a long straight circular pipe at its downstream to characterise recirculation zone and skin friction coefficient in the laminar regime. The flow is considered linearly stable for Reynolds numbers sufficiently below natural transition. A perturbation is added to the inlet fully developed laminar velocity profile to test the flow response to finite amplitude disturbances and to characterise sub-critical transition. Design/methodology/approach Direct numerical simulations of the Navier–Stokes equations have been solved using a spectral element method. Findings It is found that the onset of disordered motion and the dynamics of the localised turbulence patch are controlled by the Reynolds number, the perturbation amplitude and the half angle of the pipe. Originality/value The authors clarify different stages of flow behaviour under the finite amplitude perturbations and shed more light to flow physics such as existence of Kelvin–Helmholtz instabilities as well as mechanism of turbulent puff shedding in diverging pipe flows.

Stochastic theory and direct numerical simulations of the relative motion of high-inertia particle pairs in isotropic turbulence

2017 ◽  
Vol 813 ◽  
pp. 205-249 ◽  
Author(s):  
Rohit Dhariwal ◽  
Sarma L. Rani ◽  
Donald L. Koch
Keyword(s):  
Numerical Simulations ◽  
Relative Motion ◽  
Relative Velocity ◽  
Isotropic Turbulence ◽  
Stochastic Theory ◽  
Stokes Number ◽  
Time Correlation ◽  
Direct Numerical Simulations ◽  
Time Integral ◽  
Particle Pairs

The relative velocities and positions of monodisperse high-inertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al., J. Fluid Mech., vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian two-time correlation of fluid relative velocities seen by pairs that are nearly stationary. The two-time correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by large-scale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integral-scale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS. LS consist of evolving the Langevin equations for pair separation and relative velocity, which is statistically equivalent to solving the classical Fokker–Planck form of the pair PDF equation. Langevin simulations of particle-pair dispersion were performed using three closure forms of the diffusivity – i.e. the one containing the time integral of the Eulerian two-time correlation of the seen fluid relative velocities and the two analytical diffusivity expressions. In the first closure form, the two-time correlation was computed using DNS of forced isotropic turbulence laden with stationary particles. The two analytical closure forms have the advantage that they can be evaluated using a model for the turbulence energy spectrum that closely matched the DNS spectrum. The three diffusivities are analysed to quantify the effects of the approximations made in deriving them. Pair relative-motion statistics obtained from the three sets of Langevin simulations are compared with the results from the DNS of (moving) particle-laden forced isotropic turbulence for $St_{\unicode[STIX]{x1D702}}=10,20,40,80$ and $Re_{\unicode[STIX]{x1D706}}=76,131$. Here, $St_{\unicode[STIX]{x1D702}}$ is the particle Stokes number based on the Kolmogorov time scale and $Re_{\unicode[STIX]{x1D706}}$ is the Taylor micro-scale Reynolds number. Statistics such as the radial distribution function (RDF), the variance and kurtosis of particle-pair relative velocities and the particle collision kernel were computed using both Langevin and DNS runs, and compared. The RDFs from the stochastic runs were in good agreement with those from the DNS. Also computed were the PDFs $\unicode[STIX]{x1D6FA}(U|r)$ and $\unicode[STIX]{x1D6FA}(U_{r}|r)$ of relative velocity $U$ and of the radial component of relative velocity $U_{r}$ respectively, both PDFs conditioned on separation $r$. The first closure form, involving the Eulerian two-time correlation of fluid relative velocities, showed the best agreement with the DNS results for the PDFs.

Sign in / Sign up

Export Citation Format

Share Document