Direct simulations of the turbulent shear layer are performed for subsonic to supersonic Mach numbers. Fully developed turbulence is achieved with profiles of mean velocity and turbulence intensities that agree well with laboratory experiments. The thickness growth rate of the shear layer exhibits a large reduction with increasing values of the convective Mach number, Mc. In agreement with previous investigations, it is found that the normalized pressure–strain term decreases with increasing Mc, which leads to inhibited energy transfer from the streamwise to cross-stream fluctuations, to the reduced turbulence production observed in DNS, and, finally, to reduced turbulence levels as well as reduced growth rate of the shear layer. An analysis, based on the wave equation for pressure, with supporting DNS is performed with the result that the pressure–strain term decreases monotonically with increasing Mach number. The gradient Mach number, which is the ratio of the acoustic time scale to the flow distortion time scale, is shown explicitly by the analysis to be the key quantity that determines the reduction of the pressure–strain term in compressible shear flows. The physical explanation is that the finite speed of sound in compressible flow introduces a finite time delay in the transmission of pressure signals from one point to an adjacent point and the resultant increase in decorrelation leads to a reduction in the pressure–strain correlation.The dependence of turbulence intensities on the convective Mach number is investigated. It is found that all components decrease with increasing Mc and so does the shear stress.DNS is also used to study the effect of different free-stream densities parameterized by the density ratio, s = ρ2/ρ1, in the high-speed case. It is found that changes in the temporal growth rate of the vorticity thickness are smaller than the changes observed in momentum thickness growth rate. The momentum thickness growth rate decreases substantially with increasing departure from the reference case, s = 1. The peak value of the shear stress, uv, shows only small changes as a function of s. The dividing streamline of the shear layer is observed to move into the low-density stream. An analysis is performed to explain this shift and the consequent reduction in momentum thickness growth rate.

Large-eddy simulations of the Richtmyer-Meshkov instability with reshock are presented and the results are compared with experiments. Several configurations of shocks initially travelling from light (air) to heavy (sulfur hexafluoride, SF 6) have been simulated to match previous experiments and good agreement is found in the growth rates of the turbulent mixing zone (TMZ). The stretched-vortex subgrid model used in this study allows for subgrid continuation modelling, where statistics of the unresolved scales of the flow are estimated. In particular, this multiscale modelling allows the anisotropy of the flow to be extended to the dissipation scale, η, and estimates to be formed for the subgrid probability density function of the mixture fraction of air/SF 6 based on the subgrid variance, including the effect of Schmidt number.

No abstract

A new, simple, and computationally efficient interface capturing scheme based on a diffuse interface approach is presented for simulation of compressible multiphase flows. Multi-fluid interfaces are represented using field variables (interface functions) with associated transport equations that are augmented, with respect to an established formulation, to enforce a selected interface thickness. The resulting interface region can be set just thick enough to be resolved by the underlying mesh and numerical method, yet thin enough to provide an efficient model for dynamics of well-resolved scales. A key advance in the present method is that the interface regularization is asymptotically compatible with the thermodynamic mixture laws of the mixture model upon which it is constructed. It incorporates first-order pressure and velocity non-equilibrium effects while preserving interface conditions for equilibrium flows, even within the thin diffused mixture region. We first quantify the improved convergence of this formulation in some widely used one-dimensional configurations, then show that it enables fundamentally better simulations of bubble dynamics. Demonstrations include both a spherical bubble collapse, which is shown to maintain excellent symmetry despite the Cartesian mesh, and a jetting bubble collapse adjacent a wall. Comparisons show that without the new formulation the jet is suppressed by numerical diffusion leading to qualitatively incorrect results.

A three-dimensional direct numerical simulation (DNS) study of a spatially evolving planar turbulent reacting jet is reported. Combustion of methane with air is modelled using a four-step reduced mechanism in the non-premixed regime. A total of eight chemical species are integrated in time along with the fluid mechanical fields. The solution of the compressible Navier-Stokes equations is obtained numerically for moderately low Mach number. A large computational grid, with 100 million grid points, is required in order to resolve the flame. The cold flow Reynolds number is 3000. The focus of the study is to investigate the dynamics of extinction fronts in threedimensional turbulent flows. A novel data reduction and identification algorithm was developed to postprocess the large DNS database and extract the shape of the evolving flame surface including its edges and their propagation velocity. The joint probability density function (p.d.f.) of edge velocity and scalar dissipation was obtained and the results indicate that the three-dimensional flame edges propagate with a velocity that is largely controlled by the local rate of scalar dissipation, or equivalently in terms of the local Damköhler number at the flame edge, as predicted by theory. Naturally, the effects of unsteadiness in this flow produce a broad joint p.d.f. The statistics collected also suggest that the mean value of the hydrogen radical reaction rate are very small in the turbulent regions of the flow owing to the functional form of the hydrogen radical reaction rate itself. The consequence of these results in the context of turbulent combustion modelling is discussed. Additional statistical and morphological information of the flame is provided.