Supersonic turbulence in shock-bound interaction zones. I. Symmetric settings

Colliding hypersonic flows play a decisive role in many astrophysical objects. They contribute, for example, to the molecular cloud structure, the X-ray emission of O-stars, differentiation of galactic sheets, appearance of wind-driven structures, or, possibly, to the prompt emission of γ-ray bursts. Our intention is thorough investigation of the turbulent interaction zone of such flows, the cold dense layer (CDL). In this paper, we focus on the idealized model of a 2D plane parallel isothermal slab and on symmetric settings, where both flows have equal parameters. We performed a set of high-resolution simulations with upwind Mach-numbers, 5 < M_u < 90. We find that the CDL is irregularly shaped and has a patchy and filamentary interior. The size of these structures increases with ℓ_cdl, the extension of the CDL. On average, but not at each moment, the solution is nearly self-similar and only depends on M_u. We give the corresponding analytical expressions, with numerical constants derived from the simulation results. In particular, we find the root-mean-square Mach-number to scale as M_rms ≈ 0.2~M_u. The mean density, ρ_m ≈ 30~ρ_u is independent of M_u. The fraction f_eff of the upwind kinetic energy that survives shock passage scales as f_eff= 1 - M_rms^-0.6. This dependence persists if the upwind flow parameters differ from one side to the other of the CDL, indicating that the turbulence within the CDL and its driving are mutually coupled. Another finding points in the same direction, namely that the auto-correlation length of the confining shocks and the characteristic length scale of the turbulence within the CDL are proportional. Larger upstream Mach-numbers lead to a faster expanding CDL, confining interfaces that are less inclined with respect to the upstream flow direction, more efficient driving, and finer interior structure with respect to the extension of the CDL.