The size distribution of supernova remnants (SNRs) can help to clarify the various aspects of their evolution and interaction with the interstellar medium (ISM). Since the observed samples of SNRs are a collection of objects with very different ages and origin that evolve in different conditions of the ISM, statistical Monte Carlo methods can be used to model their statistical distributions. Based on very general assumptions on the evolution, we have modeled samples of SNRs at various initial and environmental conditions, which were then compared with observed collections of SNRs. In the evolution of SNRs the pressure of the ISM is taken into account, which determines their maximum sizes and lifetimes. When comparing the modeled and observed distributions, it is very important to have homogeneous observational data free from selection effects. We found that a recently published collection of SNRs in M33 (Long et al. 2010, ApJS,187,495) satisfies this requirement if we select the X-ray SNRs with hardness ratios in a limited range of values. An excellent agreement between distributions of this subset of SNRs and the subset of modeled SNRs was reached for a volume filling-factor of the warm phase of the ISM (partly ionized gas with $n_{\rm H}\sim 0.2-0.5~ \rm {cm}^{-3}; T \sim 8000-10000~K $) in M33 of $\sim\ 90%$. The statistical distributions constructed in this way, which reproduce practically all the statistical properties of observed SNRs, allowed us to obtain one of the important parameters of M33: the birthrate is one SNR every $ {140} - {150}$ yr, and the total number of SNRs with a shock Mach number $M_{s} \geq 2$ is larger than $\sim 1000$.

The evolution of the radio emission of shell-type Supernova remnants (SNRs) is modeled within the framework of the simple and commonly used assumptions that the mechanism of diffusive shock acceleration (DSA) is responsible for generating radio emitting electrons and that the magnetic field is the typical interstellar field compressed at the shock. It is considered that electrons are injected into the mechanism in test-particle regime directly from the high energy tail of the downstream Maxwellian distribution function. The model can be applied to most of the observed SNRs. It is shown that the model successfully explains the many averaged observational properties of evolved shell-type SNRs. In particular, the radio surface brightness ($\Sigma$) evolves with diameter as $\sim D^{-(0.3 \div 0.5)}$, while the bounding shock is strong (Mach number is ${\mathcal M} \geq10$), followed by steep decrease (steeper than $\sim D^{-4.5}$) for ${\cal M} <10$. Such evolution of the surface brightness with diameter and its strong dependence on the environmental parameters strongly reduce the usefulness of $\Sigma - D$ relations as a tool for determining the distances to SNRs. The model predicts no radio emission from SNRs in the late radiative stage of evolution and the existence of radio-quiet but relatively active SNRs is possible. Our model easily explains very large-diameter radio sources such as the Galactic Loops and the candidates for Hypernova radio remnants. The model predicts that most of the observed SNRs are located in a tenuous phase of the ISM. From the comparison of the model results with the statistics of evolved shell-type SNRs, we were able to estimate the fraction of electrons accelerated from the thermal pool in the range $(3\div 11) \times 10^{- 4}$.

Possibility of realization of scenario for the AXPs and SGRs origins from radio pulsars, which were undergone frequent and strong glitches, was analyzed. It is shown that the characteristics of such pulsars - the possible progenitors of AXPs and SGRs - their association with supernova remnants (SNRs) and evolution on the $P - \dot P$ diagram, taking into account their real ages, conflict with offered scenario.

In present study I examine the capability of diffusive shock acceleration mechanism to explain existing data on radio emission from evolved large diameter shell-type adiabatic supernova remnants (SNRs). Time-dependent ''onion-shell'' model for the radio emission of SNRs is developed, which is based on the assumptions: a) acceleration takes place from thermal energies and test-particle approximation is valid; b) the problem of injection is avoided by introducing, like Bell (1978), two injection parameters; c) to take into consideration very late stages of SNR evolution the analytic approximation of Cox and Andersen (1982) for the shell structure is used; c)no radiative cooling. Constructed Surface Brightness - Diameter $(\Sigma -D)$ tracks are compared with the empirical $\Sigma -D$ diagram. The main conclusion of the study is that the DSA mechanism is capable of explaining all the statistics of radio SNRs including very large diameter remnants and giant galactic loops.