This review aims at providing an up-to-date status and a general introduction to the subject of the numerical study of energetic particle acceleration and transport in turbulent astrophysical flows. The subject is also complemented by a short overview of recent progresses obtained in the domain of laser plasma experiments. We review the main physical processes at the heart of the production of a non-thermal distribution in both Newtonian and relativistic astrophysical flows, namely the first and second order Fermi acceleration processes. We also discuss shock drift and surfing acceleration, two processes important in the context of particle injection in shock acceleration. We analyze with some details the particle-in-cell (PIC) approach used to describe particle kinetics. We review the main results obtained with PIC simulations in the recent years concerning particle acceleration at shocks and in reconnection events. The review discusses the solution of Fokker-Planck problems with application to the study of particle acceleration at shocks but also in hot coronal plasmas surrounding compact objects. We continue by considering large scale physics. We describe recent developments in magnetohydrodynamic (MHD) simulations. We give a special emphasize on the way energetic particle dynamics can be coupled to MHD solutions either using a multi-fluid calculation or directly coupling kinetic and fluid calculations. This aspect is mandatory to investigate the acceleration of particles in the deep relativistic regimes to explain the highest Cosmic Ray energies.

Weakly magnetized, relativistic collisionless shock waves are not only the natural offsprings of relativistic jets in high-energy astrophysical sources, they are also associated with some of the most outstanding displays of energy dissipation through particle acceleration and radiation. Perhaps their most peculiar and exciting feature is that the magnetized turbulence that sustains the acceleration process, and (possibly) the secondary radiation itself, is self-excited by the accelerated particles themselves, so that the phenomenology of these shock waves hinges strongly on the microphysics of the shock. In this review, we draw a status report of this microphysics, benchmarking analytical arguments with particle-in-cell simulations, and extract consequences of direct interest to the phenomenology, regarding in particular the so-called microphysical parameters used in phenomenological studies.

Perpendicular relativistic ($\gamma_0=10$) shocks in magnetized pair plasmas are investigated using two dimensional Particle-in-Cell simulations. A systematic survey, from unmagnetized to strongly magnetized shocks, is presented accurately capturing the transition from Weibel-mediated to magnetic-reflection-shaped shocks. This transition is found to occur for upstream flow magnetizations $10^{-3}<\sigma<10^{-2}$ at which a strong perpendicular net current is observed in the precursor, driving the so-called current-filamentation instability. The global structure of the shock and shock formation time are discussed. The MHD shock jump conditions are found in good agreement with the numerical results, except for $10^{-4} < \sigma < 10^{-2}$ where a deviation up to 10\% is observed. The particle precursor length converges toward the Larmor radius of particles injected in the upstream magnetic field at intermediate magnetizations. For $\sigma>10^{-2}$, it leaves place to a purely electromagnetic precursor following from the strong emission of electromagnetic waves at the shock front. Particle acceleration is found to be efficient in weakly magnetized perpendicular shocks in agreement with previous works, and is fully suppressed for $\sigma > 10^{-2}$. Diffusive Shock Acceleration is observed only in weakly magnetized shocks, while a dominant contribution of Shock Drift Acceleration is evidenced at intermediate magnetizations. The spatial diffusion coefficients are extracted from the simulations allowing for a deeper insight into the self-consistent particle kinematics and scale with the square of the particle energy in weakly magnetized shocks. These results have implications for particle acceleration in the internal shocks of AGN jets and in the termination shocks of Pulsar Wind Nebulae.