Energy spectra of galactic cosmic rays calculated within the framework of the NoRD (nonlocal relativistic diffusion) model are presented. The model accounts for the turbulent character of the interstellar medium and the relativistic speed limit requirement. Calculations account for spiral distribution of sources, boundedness of halo, spallation of nuclei, energy dependence of the diffusion coefficient, tempered power law injection spectrum. We show that the knee in the background spectra from the ensemble of supernovae can arise due to relativistic speed limit requirement. We compare our calculations with the results obtained in frames of the local models (LoD and LoRD) and nonlocal nonrelativistic model (NoD).

Two models of anomalous diffusion of cosmic ray in the leaky-box approximation are compared: one of them is based on the decoupled time-space L\'evy flights and the other on fractional walks with a finite free motion velocity. Distributions of first passage time and paths are computed and evolution of diffusion packets to equilibrium state is shown. Calculations demonstrate essential difference between the two models: the coupled scheme gives more realistic results.

Fractional differential approach to cosmic ray physics problems is discussed. A short review in this field is given, some results are represented, analyzed and criticized. A new model called the bounded anomalous diffusion model is offered. Its equation includes the fractional material derivative which allows to take into account the finite speed of cosmic rays particles.