A comprehensive two-phase compressible gas-solid particle computational fluid dynamics-based modeling method has been developed and applied to study the erosion and fouling along flow path surfaces due to particle impact. Extending upon methods to predict particle interaction phenomena which incorporate rolling, twisting, sliding, and adhesion forces to capture conditions that frequently occur as particle laden flow passes over a flow path surface, the simulation methods have been specifically designed for high particle concentrations and large particle sizes. In the developed technique, individual particles are directly discretized in the computational mesh, and the particles move through the fluid, and interact with the flow, flow path boundaries, and other particles. With particle motion calculated from the direct two-phase flow simulations, the distribution and intensity of the mass loss or erosion conditions along a flow path surface resulting from particle impact is predicted. The developed method was applied to study the erosion of an aluminum wall due to impacting sand particles of different sizes, impact velocities, and counts. A more fluid-like particle motion and distribution upon impact were found as the particle size decreased. The method was also implemented to examine the erosion that develops for three common flow configurations which force particle-flow path surface interactions: a diverging/converging flow path, a bypass flow path, and a baffle configuration. These applications demonstrate the utility of the model to explore and better understand the relationships between the geometry, the flow, the particles, and the adhesion and erosion that develop as particle laden flow moves over the flow path surfaces.