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

An open and parallel multiresolution framework using block-based adaptive grids

A numerical approach for solving evolutionary partial differential equations in two and three space dimensions on block-based adaptive grids is presented. The numerical discretization is based on high-order, central finite-differences and explicit time integration. Grid refinement and coarsening are triggered by multiresolution analysis, i.e. thresholding of wavelet coefficients, which allow controlling the precision of the adaptive approximation of the solution with respect to uniform grid computations. The implementation of the scheme is fully parallel using MPI with a hybrid data structure. Load balancing relies on space filling curves techniques. Validation tests for 2D advection equations allow to assess the precision and performance of the developed code. Computations of the compressible Navier-Stokes equations for a temporally developing 2D mixing layer illustrate the properties of the code for nonlinear multi-scale problems. The code is open source

Self-Similar Evolution of Cosmic-Ray-Modified Quasi-Parallel Plane Shocks

Using an improved version of the previously introduced CRASH (Cosmic Ray Acceleration SHock) code, we have calculated the time evolution of cosmic-ray (CR) modified quasi-parallel plane shocks for Bohm-like diffusion, including self-consistent models of Alfven wave drift and dissipation, along with thermal leakage injection of CRs. The new simulations follow evolution of the CR distribution to much higher energies than our previous study, providing a better examination of evolutionary and asymptotic behaviors. The postshock CR pressure becomes constant after quick initial adjustment, since the evolution of the CR partial pressure expressed in terms of a momentum similarity variable is self-similar. The shock precursor, which scales as the diffusion length of the highest energy CRs, subsequently broadens approximately linearly with time, independent of diffusion model, so long as CRs continue to be accelerated to ever-higher energies. This means the nonlinear shock structure can be described approximately in terms of the similarity variable, x/(u_s t), where u_s is the shock speed once the postshock pressure reaches an approximate time asymptotic state. As before, the shock Mach number is the key parameter determining the evolution and the CR acceleration efficiency, although finite Alfven wave drift and wave energy dissipation in the shock precursor reduce the effective velocity change experienced by CRs, so reduce acceleration efficiency noticeably, thus, providing a second important parameter at low and moderate Mach numbers.Comment: 29 pages, 8 figure

Efficiency of Nonlinear Particle Acceleration at Cosmic Structure Shocks

We have calculated the evolution of cosmic ray (CR) modified astrophysical shocks for a wide range of shock Mach numbers and shock speeds through numerical simulations of diffusive shock acceleration (DSA) in 1D quasi- parallel plane shocks. The simulations include thermal leakage injection of seed CRs, as well as pre-existing, upstream CR populations. Bohm-like diffusion is assumed. We model shocks similar to those expected around cosmic structure pancakes as well as other accretion shocks driven by flows with upstream gas temperatures in the range $T_0=10^4-10^{7.6}$K and shock Mach numbers spanning $M_s=2.4-133$. We show that CR modified shocks evolve to time-asymptotic states by the time injected particles are accelerated to moderately relativistic energies (p/mc \gsim 1), and that two shocks with the same Mach number, but with different shock speeds, evolve qualitatively similarly when the results are presented in terms of a characteristic diffusion length and diffusion time. For these models the time asymptotic value for the CR acceleration efficiency is controlled mainly by shock Mach number. The modeled high Mach number shocks all evolve towards efficiencies $\sim 50$%, regardless of the upstream CR pressure. On the other hand, the upstream CR pressure increases the overall CR energy in moderate strength shocks ($M_s \sim {\rm a few}$). (abridged)Comment: 23 pages, 12 ps figures, accepted for Astrophysical Journal (Feb. 10, 2005

Study of the microstructure, tensile properties and hardness of AZ61 magnesium alloy subjected to severe plastic deformation

Hot extruded (EX) AZ61 magnesium alloy was processed by the twist channel angular pressing (TCAP) method, which combines equal channel angular pressing (ECAP) and twist extrusion (TE) processes and significantly improves the efficiency of the grain refinement process. Both the initial hot extruded AZ61 alloy and the alloy after completion of TCAP processing were examined by using optical microscopy (OM), scanning electron microscopy (SEM), transmission electron microscopy (TEM) and electron backscatter diffraction (EBSD) and their corresponding micro-tensile testing (M-TT) and hardness testing at room temperature. The results showed that the microstructure of hot extruded alloy was refined well by TCAP due to dynamic recrystallization (DRX) caused by TCAP. The tensile properties, investigated by micro-tensile testing (M-TT), of the AZ61 alloy were significantly improved due to refined microstructure. The highest tensile properties including YS of 240.8 MPa, UTS of 343.6 MPa and elongation of 21.4% of the fine-grained alloy with average grain size below 1.5 mu m was obtained after the third TCAP pass at 200 degrees C using the processing route B-c.Web of Science810art. no. 77