The development of a parallel adaptive mesh refinement (AMR) scheme is described for solving the governing equations for multi-phase (gas–particle) core flows in solid propellant rocket motors (SRMs). An Eulerian formulation is used to describe the coupled motion between the gas and particle phases. A cell-centred upwind finite-volume discretization and the use of limited linear reconstruction, Riemann solver based flux functions for the gas and particle phases, and explicit multi-stage time-stepping allows for high solution accuracy and computational robustness. A Riemann problem is formulated for prescribing boundary data at the burning surface and a mesh adjustment algorithm has been implemented to adjust the multi-block quadrilateral mesh to the combustion interface. A flexible block-based hierarchical data structure is used to facilitate automatic solution-directed mesh adaptation according to physics-based refinement criteria. Efficient and scalable parallel implementations are achieved with domain decomposition on distributed memory multi-processor architectures. Numerical results are described to demonstrate the capabilities of the approach for predicting SRM core flows.