The modeling assumptions and governing equations required to estimate the heatshield surface recession due to dust particle impacts are presented. The dust particle size distribution can be estimated using a modified gamma distribution. The vertical distribution of dust particles in the atmosphere in terms of the ratio of dust particle to atmospheric density can be related to measured or computed values of the dust opacity. The different approaches to coupling the particle and fluid dynamics are discussed. The equations and modeling assumptions presented in this paper are applied to computing the surface recession due to dust particle impacts on the ExoMars Schiaparelli capsule entering the Martian atmosphere during a representative dust storm. The dust surface recession at the stagnation point of about 1 millimeter is about 40 percent of the recession value due to charring ablation.