We investigate the shock-induced turbulent mixing between a light and a heavy gas, where a Richtmyer-Meshkov instability (RMI) is initiated by a shock wave with Mach number Ma = 1.5. The prescribed initial conditions define a deterministic multimode interface perturbation between the gases, which can be imposed exactly for different simulation codes and resolutions to allow for quantitative comparison. Wellresolved large-eddy simulations are performed using two different and independently developed numerical methods with the objective of assessing turbulence structures, prediction uncertainties and convergence behaviour. The two numerical methods differ fundamentally with respect to the employed subgrid-scale regularisation, each representing state-of-the-art approaches to RMI. Unlike previous studies, the focus of the present investigation is to quantify the uncertainties introduced by the numerical method, as there is strong evidence that subgrid-scale regularisation and truncation errors may have a significant effect on the linear and nonlinear stages of the RMI evolution. Fourier diagnostics reveal that the larger energy-containing scales converge rapidly with increasing mesh resolution and thus are in excellent agreement for the two numerical methods. Spectra of gradient-dependent quantities, such as enstrophy and scalar dissipation rate, show stronger dependences on the small-scale flow field structures as a consequence of truncation error effects, which for one numerical method are dominantly dissipative and for the other dominantly dispersive. Additionally, the study reveals details of various stages of RMI, as the flow transitions from large-scale nonlinear entrainment to fully developed turbulent mixing. The growth rates of the mixing zone widths as obtained by the two numerical methods are ~t^7/12 before re-shock and ~(t - t₀)^2/7 long after re-shock. The decay rate of turbulence kinetic energy is consistently ~(t-t₀)^-10/7 at late times, where the molecular mixing fraction approaches an asymptotic limit ≈ 0.85. The anisotropy measure.

We quantify initial-data uncertainties on a shock accelerated heavy-gas cylinder by two-dimensional well-resolved direct numerical simulations. A high-resolution compressible multicomponent flow simulation model is coupled with a polynomial chaos expansion to propagate the initial-data uncertainties to the output quantities of interest. The initial flow configuration follows previous experimental and numerical works of the shock accelerated heavy-gas cylinder.We investigate three main initialdata uncertainties, (i) shock Mach number, (ii) contamination of SF₆ with acetone, and (iii) initial deviations of the heavy-gas region from a perfect cylindrical shape. The impact of initial-data uncertainties on the mixing process is examined. The results suggest that the mixing process is highly sensitive to input variations of shock Mach number and acetone contamination. Additionally, our results indicate that the measured shock Mach number in the experiment of Tomkins et al. ["An experimental investigation of mixing mechanisms in shock-accelerated flow," J. Fluid. Mech. 611, 131 (2008)] and the estimated contamination of the SF₆ region with acetone [S. K. Shankar, S. Kawai, and S. K. Lele, "Two-dimensional viscous flow simulation of a shock accelerated heavy gas cylinder," Phys. Fluids 23, 024102 (2011)] exhibit deviations from those that lead to best agreement between our simulations and the experiment in terms of overall flow evolution.

We present results of well-resolved direct numerical simulations (DNS) of the turbulent flow evolving from Richtmyer-Meshkov instability (RMI) in a shock-tube with square cross section. The RMI occurs at the interface between a mixture of O2 and N2 (light gas) and SF6 and acetone (heavy gas). The interface between the light and heavy gas is accelerated by a Ma = 1.5 planar shock wave. RMI is triggered by a well-defined multimodal initial disturbance at the interface. The DNS exhibit grid-resolution independent statistical quantities and support the existence of a Kolmogorov-like inertial range with a k-5/3 scaling unlike previous simulations found in the literature. The results are in excellent agreement with the experimental data of Weber et al. ["Turbulent mixing measurements in the Richtmyer-Meshkov instability," Phys. Fluids24, 074105 (2012)]10.1063/1.4733447.

In most technical applications involving cavitation, vapor bubbles occur in clouds, and their collapse is affected by the interaction with neighboring bubbles. One approach to study the influence of these interactions is the investigation of the collapse of cavity arrays in water under shock wave loading. We describe in detail the collapse mechanisms during the collapse of a horizontal cavity array, with particular consideration of maximum pressures. As general trend, we find a pressure amplification in consecutive cavity collapses. However, by increasing the number of cavities, we are able to demonstrate that the amplification is not monotonic. A parameter study of the bubble separation distance in horizontal arrays shows that a smaller distance generally, but not necessarily, results in larger collapse pressure. Exceptions from the general trend are due to the very complex shock and expansion-wave interactions and demonstrate the importance of using state-of-the-art numerical methods. By varying boundary conditions, we illustrate the significance of large test sections in experimental investigations, as the expansion wave emitted at a free surface has a large effect on the collapse dynamics.

In this paper, we investigate the high-speed dynamics of symmetric and asymmetric cavitation bubble-collapse. For this purpose, a sharp-interface numerical model is employed, that includes a numerically efficient evaporation/condensation model. The underlying assumption is that phase change occurs in thermal non-equilibrium and that the associated timescale is much larger than that of the wave-dynamics described by the interfacial Riemann problem. The sharp-interface model allows for an accurate tracking of the interface evolution throughout collapse and rebound. With a first set of simulations, we investigate the influence of the non-equilibrium on the relaxation behaviour of an oscillating vapour bubble. We observe that a good prediction of the phase-change rate is essential. Of high practical interest is the collapse of cavitation bubbles near walls under high ambient-pressure conditions. We investigate the differences in collapse evolution for detached and attached bubbles. It is shown that the maximum wall pressure strongly depends on the symmetry of the collapse mechanisms, and regions with a high probability of bubble rebound are identified. Asymmetric attached bubbles lead to significantly different topology changes during collapse than symmetric bubbles but exhibit roughly the same range of maximum pressures.

We propose a conservative, second-order accurate immersed interface method for representing incompressible fluid flows over complex three dimensional solid obstacles on a staggered Cartesian grid. The method is based on a finite-volume discretization of the incompressible Navier-Stokes equations which is modified locally in cells that are cut by the interface in such a way that accuracy and conservativity are maintained. A level-set technique is used for description and tracking of the interface geometry, so that an extension of the method to moving boundaries and flexible walls is straightforward. Numerical stability is ensured for small cells by a conservative mixing procedure. Discrete conservation and sharp representation of the fluid-solid interface render the method particularly suitable for Large-Eddy Simulations of high-Reynolds number flows. Accuracy, second-order grid convergence and robustness of the method is demonstrated for several test cases: inclined channel flow at Re= 20, flow over a square cylinder at Re= 100, flow over a circular cylinder at Re= 40, Re= 100 and Re= 3900, as well as turbulent channel flow with periodic constrictions at Re= 10,595.