A computational method of simulations for processes of heterogeneous hydrodynamics with take of phase transitions will be discussed. The method is based on relevant approximation of conservation laws for mass, momentum, and energy in integral and differential forms. The time and spatial approximation is natural and numerical simulations are realized as direct computer experiments. It is supposed that the fluids are compressible and non-viscous. Heterogeneities of the fluids are considered as small drops or particles of one fluid within other fluid. Total number of the drops may be large enough and the drops may have phase transitions. Therefore, simulations of the main fluid with small transited drops dynamics are considered. The particle dynamics will be modelled as in the particle-in-cell method, and in the main fluid as in the large particle method. This approach makes it possible to simulate phase transitions under certain assumptions about heterogeneous fluids. The calculation algorithm of this method is implemented as a computer simulation of the dynamics of a multiphase carrier fluid containing particles that can undergo, for example, graphite-diamond phase transitions. Such transitions are modelled on the basis of the theory of phase transformations and the laws of thermodynamics. In fact, the method is a combination of the Harlow's particle-in-cell method, Belotserkovskii's large particles method and Bakhvalov's homogenization method. A modification of this method has also been developed to take into account the effects of viscosity when simulating the dynamics of a multiphase fluid in porous media. A model of the motion of such a liquid in a porous medium is obtained by freezing the motion of particles of the corresponding size in the presented method. The method will certainly be promising for numerical simulations of absorption and diffusion processes in complex fluids with phase transitions.