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.

Possible formulations of the problem of cosmic rays acceleration in the interstellar galactic medium are considered with the use of fractional differential equations. The applied technique has been physically justified. A Fermi result has been generalized to the case of the acceleration of particles in shock waves in the supernovae remnants fractally distributed in the Galaxy.

We consider the propagation of galactic cosmic rays under assumption that the interstellar medium is a fractal one. An anomalous diffusion equation in terms of fractional derivatives is used to describe of cosmic ray propagation. The anomaly in used model results from large free paths ("Levy flights") of particles between galactic inhomogeneities and long time a particle stays in a trap. An asymptotical solution of the anomalous diffusion equation for point instantaneous and impulse sources with inverse power spectrum relating to supernova bursts is found. It covers both subdiffusive and superdiffusive regimes and is expressed in terms of the stable distributions. The energy dependence of spectral exponent of observed particles in different regimes is discussed.