The aim of this work is to construct an algorithm for visualizing a polydisperse phase of solid particles (dust) in an inhomogeneous flow of a two-phase gas-dust mixture that would allow us to see, within one plot, the degree of polydispersity of the dust phase and the difference in the spatial distributions of individual fractions of dust particles in the computational domain. The developed technique allows us to reproduce concentrations from one to three fractions of dust particles in each cell in the computational domain. Each of the three fractions of dust particles is mapped to one of the main channels of the RGB palette. The intensity of the color shade is set to be proportional to the relative concentration of dust particles in this fraction. The final image for a polydisperse mixture is obtained by adding images in each of the three color channels. To visualize the degree of polydispersity, I propose depicting the spatial distribution of the entropy of the dust mixture. The definition of the entropy of a mixture is generalized to take into account the states of a mixture with zero number of particles in the mixture. They correspond to dust-free sections of the computational domain (voids). The proposed method for visualizing the polydispersity of a mixture of particles is demonstrated using the example of dynamic numerical modeling of the spatial features of dust structures formed in turbulent gas-dust flows and in flows with shock waves.

We discuss an idea whether spherical blast waves can amplify by a non-local resonant hydrodynamic mechanism inhomogeneities formed by turbulence or phase segregation in the interstellar medium. We consider the problem of a blast-wave-turbulence interaction in the Linear Interaction Approximation. Mathematically, this is an eigenvalue problem for finding the structure and amplitude of eigenfunctions describing the response of the shock-wave flow to forced oscillations by external perturbations in the ambient interstellar medium. Linear analysis shows that the blast wave can amplify density and vorticity perturbations for a wide range of length scales with amplification coefficients of up to 20, with amplification the greater, the larger the length. There also exist resonant harmonics for which the gain becomes formally infinite in the linear approximation. Their orbital wavenumbers are within the range of macro- ($l \sim 1$), meso- ($l \sim 20$) and microscopic ($l > 200$) scales. Since the resonance width is narrow: typically, $\Delta l <1$, resonance should select and amplify discrete isolated harmonics. We speculate as to a possible explanation of an observed regular filamentary structure of regular-shaped round supernova remnants such as SNR 1572, 1006 or 0509-67.5. Resonant mesoscales found ($l \approx 18 $) are surprisingly close to the observed scales ($l \approx 15$) of ripples in the shell's surface of SNR 0509-67.5.

The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO interpolation. The time integration is carried out using the third order TVD method of Runge-Kutta type, thus ensuring overall third order accuracy on smooth solutions. The magnetic field is kept near divergence-free by means of the method of generalized Lagrange multiplier. The test simulations, which include linear and non-linear continuous plasma waves, shock waves, strong explosions and the tearing instability, show that the scheme is sufficiently robust and confirm its accuracy.