Simulations of cosmic ray propagation.

We review numerical methods for simulations of cosmic ray (CR) propagation on galactic and larger scales. We present the development of algorithms designed for phenomenological and self-consistent models of CR propagation in kinetic description based on numerical solutions of the Fokker-Planck equation. The phenomenological models assume a stationary structure of the galactic interstellar medium and incorporate diffusion of particles in physical and momentum space together with advection, spallation, production of secondaries and various radiation mechanisms.

Resonant Kelvin-Helmholtz modes in sheared relativistic flows.

Certain aspects of the (linear and nonlinear) stability of sheared relativistic (slab) jets are analyzed. The linear problem has been solved for a wide range of jet models well inside the ultrarelativistic domain (flow Lorentz factors up to 20, specific internal energies approximately 60c2). As a distinct feature of our work, we have combined the analytical linear approach with high-resolution relativistic hydrodynamical simulations, which has allowed us (i) to identify, in the linear regime, resonant modes specific to the relativistic shear layer, (ii) to confirm the result of the linear analysis with numerical simulations, and (iii) more interestingly, to follow the instability development through the nonlinear regime.