The forward-reverse shock pair at large heliocentric distances

An unsteady one-dimensional numerical magnetohydrodynamic (MHD) model is developed in order to study the essential physical processes involved in the development of the forward-reverse shock pair in the heliosphere. In the model, MHD shocks are treated as boundary surfaces which divide the domain of interest in the r-t plane into several flow regions. The positions of the shock boundary surfaces between two neighboring flow regions are determined by shock speed. On the basis of integrations of the model, it is found that the strong MHD disturbances generated in a corotating interaction region (CIR) propagate at a fast speed relative to the moving material, and that the wave propagation speed is greater in CIR than in its surroundings. This causes disturbances in CIR to pile up and form a shock pair. The newly formed shock pair will in turn propagate outward from the leading edge to interact with ambient rarefaction regions. This interaction accounts for the double sawtooth configuration observed in velocity profiles of shock pairs. It is also demonstrated that the merging of two shocks produces a stronger shock and constant surface on its backside. Computer generated velocity profiles based on the model are presented.