PIC simulation of a shock tube: Implications for wave transmission in the heliospheric boundary region

A shock tube problem is solved numerically by using one-dimensional full particle-in-cell simulations under the condition that a relatively tenuous and weakly magnetized plasma is continuously pushed by a relatively dense and strongly magnetized plasma having supersonic relative velocity. A forward and a reverse shock and a contact discontinuity are self-consistently reproduced. The spatial width of the contact discontinuity increases as the angle between the discontinuity normal and ambient magnetic field decreases. The inner structure of the discontinuity shows different profiles between magnetic field and plasma density, or pressure, which is caused by a non-MHD effect of the local plasma. The region between the two shocks is turbulent. The fluctuations in the relatively dense plasma are compressible and propagating away from the contact discontinuity, although the fluctuations in the relatively tenuous plasma contain both compressible and incompressible components. The source of the compressible fluctuations in the relatively dense plasma is in the relatively tenuous plasma. Only compressible fast mode fluctuations generated in the relatively tenuous plasma are transmitted through the contact discontinuity and propagate in the relatively dense plasma. These fast mode fluctuations are steepened when passing the contact discontinuity. This wave steepening and probably other effects may cause the broadening of the wave spectrum in the very local interstellar medium plasma. The results are discussed in the context of the heliospheric boundary region or heliopause.
Comment: 14 figures, accepted for publication in ApJ