Small-amplitude perturbations in the three-dimensional cylindrical Richtmyer–Meshkov instability

We first study the linear stability of an interface between two fluids following the passage of an imploding or exploding shock wave. Assuming incompressible flow between the refracted waves following shock impact, we derive an expression for the asymptotic growth rate for a three-dimensional combination of azimuthal and axial perturbations as a function of the Atwood ratio, the axial and azimuthal wave numbers, the initial radial position and perturbation amplitude of the interface, and the interface velocity gain due to the shock interaction. From the linearized theory, a unified expression for the impulsive asymptotic growth rate in plane, cylindrical, and spherical geometries is obtained which clearly delineates the effects of perturbation growth due to both geometry and baroclinic vorticity deposition. Several different limit cases are investigated, allowing recovery of Mikaelian's purely azimuthal theory and Richtmyer's plane model. We discuss the existence of three-dimensional perturbations with zero growth, typical of curvilinear geometries, as first observed by Mikaelian. The effect of shock proximity on the interface growth rate is studied in the case of a reflected shock. Analytical predictions of the effect of the incident shock strength and the perturbation wave numbers are then compared with results obtained from highly resolved numerical simulations of cylindrical imploding Richtmyer–Meshkov instability for ideal gases. A parallel is made with the instability growth in spherical and plane geometry. In particular, we propose a representation of the perturbation growth by considering the volume of the perturbed layer. This volume is found to grow faster in the plane case than in the imploding cylindrical geometry, among other results

Atwood ratio dependence of Richtmyer-Meshkov flows under reshock conditions using large-eddy simulations

We study the shock-driven turbulent mixing that occurs when a perturbed planar density interface is impacted by a planar shock wave of moderate strength and subsequently reshocked. The present work is a systematic study of the influence of the relative molecular weights of the gases in the form of the initial Atwood ratio A. We investigate the cases A = ± 0.21, ±0.67 and ±0.87 that correspond to the realistic gas combinations air–CO_2, air–SF_6 and H_2–air. A canonical, three-dimensional numerical experiment, using the large-eddy simulation technique with an explicit subgrid model, reproduces the interaction within a shock tube with an endwall where the incident shock Mach number is ~1.5 and the initial interface perturbation has a fixed dominant wavelength and a fixed amplitude-to-wavelength ratio ~0.1. For positive Atwood configurations, the reshock is followed by secondary waves in the form of alternate expansion and compression waves travelling between the endwall and the mixing zone. These reverberations are shown to intensify turbulent kinetic energy and dissipation across the mixing zone. In contrast, negative Atwood number configurations produce multiple secondary reshocks following the primary reshock, and their effect on the mixing region is less pronounced. As the magnitude of A is increased, the mixing zone tends to evolve less symmetrically. The mixing zone growth rate following the primary reshock approaches a linear evolution prior to the secondary wave interactions. When considering the full range of examined Atwood numbers, measurements of this growth rate do not agree well with predictions of existing analytic reshock models such as the model by Mikaelian (Physica D, vol. 36, 1989, p. 343). Accordingly, we propose an empirical formula and also a semi-analytical, impulsive model based on a diffuse-interface approach to describe the A-dependence of the post-reshock growth rate

New Epicenters for Production Development in Port Cities: The Digital Innovation Hub in Genoa

In the framework of infrastructural upgrading that the port city of Genoa has been going through for at least two decades, the episode of the Erzelli Science and Technology Park represents a unicum for geographic location, functional programme, implementation process, and actors involved. Located on the hill of the eponymous name, the Park hosts the Liguria Digital Innovation Hub, responding to a need for delocalisation and territorial aggregation of large activities related to technology, production, the service sector and scientific research. The contribution explores how the realization of the Park addresses critical issues related to accessibility and to the attractiveness of the territories, declining the theme of development epicenters from a technological, productive and tertiary point of view

Turbulent mixing driven by spherical implosions. Part 1. Flow description and mixing-layer growth

We present large-eddy simulations (LES) of turbulent mixing at a perturbed, spherical interface separating two fluids of differing densities and subsequently impacted by a spherically imploding shock wave. This paper focuses on the differences between two fundamental configurations, keeping fixed the initial shock Mach number ≈1.2, the density ratio (precisely |A_0|≈0.67) and the perturbation shape (dominant spherical wavenumber ℓ_0=40 and amplitude-to-initial radius of 3%): the incident shock travels from the lighter fluid to the heavy fluid or, inversely, from the heavy to the light fluid. After describing the computational problem we present results on the radially symmetric flow, the mean flow, and the growth of the mixing layer. Turbulent statistics are developed in Part 2 (Lombardini, M., Pullin, D. I. & Meiron, D. I. J. Fluid Mech., vol. 748, 2014, pp. 113–142). A wave-diagram analysis of the radially symmetric flow highlights that the light–heavy mixing layer is processed by consecutive reshocks, and not by reverberating rarefaction waves as is usually observed in planar geometry. Less surprisingly, reshocks process the heavy–light mixing layer as in the planar case. In both configurations, the incident imploding shock and the reshocks induce Richtmyer–Meshkov (RM) instabilities at the density layer. However, we observe differences in the mixing-layer growth because the RM instability occurrences, Rayleigh–Taylor (RT) unstable scenarios (due to the radially accelerated motion of the layer) and phase inversion events are different. A small-amplitude stability analysis along the lines of Bell (Los Alamos Scientific Laboratory Report, LA-1321, 1951) and Plesset (J. Appl. Phys., vol. 25, 1954, pp. 96–98) helps quantify the effects of the mean flow on the mixing-layer growth by decoupling the effects of RT/RM instabilities from Bell–Plesset effects associated with geometric convergence and compressibility for arbitrary convergence ratios. The analysis indicates that baroclinic instabilities are the dominant effect, considering the low convergence ratio (≈2) and rather high (ℓ>10) mode numbers considered

Turbulent mixing driven by spherical implosions. Part 2. Turbulence statistics

We present large-eddy simulations (LES) of turbulent mixing at a perturbed, spherical interface separating two fluids of differing densities and subsequently impacted by a spherically imploding shock wave. This paper focuses on the differences between two fundamental configurations, keeping fixed the initial shock Mach number ≈ 1.2, the density ratio (precisely |A_0| ≈ 0.67) and the perturbation shape (dominant spherical wavenumber ℓ_0=40 and amplitude-to-initial radius of 3 %): the incident shock travels from the lighter fluid to the heavy one, or inversely, from the heavy to the light fluid. In Part 1 (Lombardini, M., Pullin, D. I. & Meiron, D. I., J. Fluid Mech., vol. 748, 2014, pp. 85-112), we described the computational problem and presented results on the radially symmetric flow, the mean flow, and the growth of the mixing layer. In particular, it was shown that both configurations reach similar convergence ratios ≈2. Here, turbulent mixing is studied through various turbulence statistics. The mixing activity is first measured through two mixing parameters, the mixing fraction parameter Theta and the effective Atwood ratio A(e), which reach similar late time values in both light-heavy and heavy-light configurations. The Taylor-scale Reynolds numbers attained at late times are estimated ≈2000 in the light-heavy case and 1000 in the heavy-light case. An analysis of the density self-correlation b, a fundamental quantity in the study of variable-density turbulence, shows asymmetries in the mixing layer and non-Boussinesq effects generally observed in high-Reynolds-number Rayleigh-Taylor (RT) turbulence. These traits are more pronounced in the light-heavy mixing layer, as a result of its flow history, in particular because of RT-unstable phases (see Part 1). Another measure distinguishing light-heavy from heavy-light mixing is the velocity-to-scalar Taylor microscales ratio. In particular, at late times, larger values of this ratio are reported in the heavy-light case. The late-time mixing displays the traits some of the traits of the decaying turbulence observed in planar Richtmyer-Meshkov (RM) flows. Only partial isotropization of the flow (in the sense of turbulent kinetic energy (TKE) and dissipation) is observed at late times, the Reynolds normal stresses (and, thus, the directional Taylor microscales) being anisotropic while the directional Kolmogorov microscales approach isotropy. A spectral analysis is developed for the general study of statistically isotropic turbulent fields on a spherical surface, and applied to the present flow. The resulting angular power spectra show the development of an inertial subrange approaching a Kolmogorov-like -5/3 power law at high wavenumbers, similarly to the scaling obtained in planar geometry. It confirms the findings of Thomas & Kares (Phys. Rev. Lett., vol. 109, 2012, 075004) at higher convergence ratios and indicates that the turbulent scales do not seem to feel the effect of the spherical mixing-layer curvature