The passage of planar shocks in a dusty gas was investigated to note effects due to particle loading and initial shock Mach number. Two-phase flow equations have been added to a conservative, monotonic flow solver to allow study of compressible particle and droplet flows, which are of importance for shock propagation in two-phase flows and spray propulsion systems. The formulation developed herein employed a conservative Eulerian treatment for the gas and particle phases. The computations were performed using the finite element method-flux corrected transport (FEM-FCT) scheme, which has shown excellent predictive capability of various compressible flows which include both strong and weak shocks. The flux limiting technique was modified to provide monotonic particle velocity fields to increase the scheme's computational stability. Adaptive unstructured methodology based on adapting to high gradients of both the fluid and particle densities was used in conjunction with the conservative shock-capturing scheme to adequately resolve strong flowfield gradients. The shock attenuation of this scheme was then compared with previous experimental and numerical results and was found to yield robust predictions. Various interphase coupling terms were also considered to note their effect on the shock attenuation.