A physical model was proposed for predicting erosion in tube and shell heat exchangers. The proposed model includes all variables and parameters that influence erosion into simulation. Erosions at tube sheet, tube ends and tubes have different relations to particle size. Particles were classified into three groups based on the dependence of erosion rate on particle size.