The current study documents the development of a new multiphase solver for the simulation of dilute dusty gas flows in the OpenFOAM suite. Even though various solvers for simulation of multiphase flows are developed in OpenFOAM, the lack of a density-based Riemann solver (which is the desirable model for the description of the interaction of shock waves with a dispersed particulate phase) is evident. As the code is developed in this opensource framework more users with broader backgrounds can use, assess and extend the code to include new features for studying new complex phenomena that are present in multiphase gas-solid flows. Moreover, the conventional dusty gas model is modified in such a way that the eigenvalue degeneracy is circumvented. The idea is that rather than devising new flux functions or modifying all the previously well-established numerical flux functions to be applicable to the pressureless Euler equation, to modify the equation model in a way that fluxes can be applied to our system. The simulation results for problems in one, two, and three-dimensional space indicate that the developed solver is a promising tool in the study of these types of flows. Furthermore, the computationally demanding problem of shock-bubble interaction in three-dimensional space is investigated with a twist. Here the bubble is replaced with a spherical dust cloud for the first time and with the help of highperformance computing facilities, justifications on topology deformation of the cloud in interaction with dust particles of different diameters and loadings are discussed in detail.