The $$YZ\beta $$YZβ shock-capturing technique, which is residual-based, was introduced in conjunction with the Streamline-Upwind/Petrov---Galerkin (SUPG) formulation of compressible flows in conservation variables. It was later also combined with the variable subgrid scale (V-SGS) formulation of compressible flows in conservation variables and successfully tested on 2D and 3D computation of inviscid flows with shocks. In this paper we extend that combined method to inviscid flow computations with particle tracking and particle---shock interaction. Particles are tracked individually, assuming one-way dependence between the particle dynamics and the flow. We present two steady-state test computations with particle---shock interaction, one in 2D and one in 3D, and show that the overall method is effective in particle tracking and particle---shock interaction analysis in compressible flows.