We examine with particle-in-cell (PIC) simulations how a parallel shock in pair plasma reacts to upstream waves, which are driven by escaping downstream particles. Initially, the shock is sustained in the two-dimensional simulation by a magnetic filamentation (beam-Weibel) instability. Escaping particles drive an electrostatic beam instability upstream. Modifications of the upstream plasma by these waves hardly affect the shock. In time, a decreasing density and increasing temperature of the escaping particles quench the beam instability. A larger thermal energy along than perpendicular to the magnetic field destabilizes the pair-Alfv\'en mode. In the rest frame of the upstream plasma, the group velocity of the growing pair-Alfv\'en waves is below that of the shock and the latter catches up with the waves. Accumulating pair-Alfv\'en waves gradually change the shock in the two-dimensional simulation from a Weibel-type shock into an Alfv\'enic shock with a Mach number that is about 6 for our initial conditions.

The problem of accelerating cosmic rays is one of fundamental importance, particularly given the uncertainty in the conditions inside the acceleration sites. Here we examine Diffusive Shock Acceleration in arbitrary turbulent magnetic fields, constructing a new model that is capable of bridging the gap between the very weak ($\delta B/B_0 \ll 1$) and the strong turbulence regimes. To describe the diffusion we provide quantitative analytical description of the "Bohm exponent" in each regime. We show that our results converge to the well known quasi-linear theory in the weak turbulence regime. In the strong regime, we quantify the limitations of the Bohm-type models. Furthermore, our results account for the anomalous diffusive behaviour which has been noted previously. Finally, we discuss the implications of our model in the study of possible acceleration sites in different astronomical objects.