We aim to elucidate cosmic ray (CR) propagation in the weakly ionized environments of supernova remnants (SNRs) basing our analysis on the cosmic ray cloud (CRC) model developed by \citet{2013ApJ...768...73M} and \citet{2016MNRAS.461.3552N}. We solved two transport equations simultaneously: one for the CR pressure and one for the Alfv\'en wave energy density where CRs are initially confined in the SNR shock. Cosmic rays trigger a streaming instability and produce slab-type resonant Alfv\'en modes. The self-generated turbulence is damped by ion-neutral collisions and by noncorrelated interaction with Alfv\'en modes generated at large scales. We show that CRs leaking in cold dense phases such as those found in cold neutral medium (CNM) and diffuse molecular medium (DiM) can still be confined over distances of a few tens of parsecs from the CRC center for a few thousand years. At 10 TeV, CR diffusion can be suppressed by two or three orders of magnitude. This effect results from a reduced ion-neutral collision damping in the decoupled regime. We calculate the grammage of CRs in these environments. We find that in both single and multi-phase setups at 10 GeV, CNM and DiM media can produce grammage in the range 10-20 $\rm{g/cm^2}$ in the CNM and DiM phases. At 10 TeV, because of nonlinear propagation the grammage increases to values in the range 0.5-20 $\rm{g/cm^2}$ in these two phases. We also present preliminary calculations in inhomogeneous interstellar medium combining two or three different phases where we obtain the same trends.

Cosmic rays (CRs) are thought to play a dynamically important role in several key aspects of galaxy evolution, including the structure of the interstellar medium, the formation of galactic winds, and the non-thermal pressure support of halos. We introduce a numerical model solving for the CR streaming instability and acceleration of CRs at shocks with a fluid approach in the adaptive mesh refinement code ramses. CR streaming is solved with a diffusion approach and its anisotropic nature is naturally captured. We introduce a shock finder for the ramses code that automatically detects shock discontinuities in the flow. Shocks are the loci for CR injection, and their efficiency of CR acceleration is made dependent on the upstream magnetic obliquity according to the diffuse shock acceleration mechanism. We show that the shock finder accurately captures shock locations and estimates the shock Mach number for several problems. The obliquity-dependent injection of CRs in the Sedov solution leads to situations where the supernova bubble exhibits large polar caps (homogeneous background magnetic field), or a patchy structure of the CR distribution (inhomogeneous background magnetic field). Finally, we combine both accelerated CRs with streaming in a simple turbulent interstellar medium box, and show that the presence of CRs significantly modifies the structure of the gas.

We study the problem of the escape and transport of Cosmic-Rays (CR) from a source embedded in a fully ionised, hot phase of the interstellar medium (HIM). In particular, we model the CR escape and their propagation in the source vicinity taking into account excitation of Alfv\'enic turbulence by CR streaming and mechanisms damping the self-excited turbulence itself. Our estimates of escape radii and times result in large values (100 pc, $2\times10^5$ yr) for particle energies $\lesssim20$ GeV and smaller values for particles with increasing energies (35 pc and 14 kyr at 1 TeV). These escape times and radii, when used as initial conditions for the CR propagation outside the source, result in relevant suppression of the diffusion coefficient (by a factor 5-10) on time-scales comparable with their (energy dependent) escape time-scale. The damping mechanisms are fast enough that even on shorter time scales the Alfv\'enic turbulence is efficiently damped, and the ratio between random and ordered component of the magnetic field is $\delta B/B_0\ll 1$, justifying the use of quasi-linear theory. In spite of the suppressed diffusion coefficient, and then the increased residence time in the vicinity (<200 pc) of their source, the grammage accumulated by CRs after their escape is found to be negligible (at all energies) as compared to the one accumulated while diffusing in the whole Galaxy, due to the low density of the HIM.