The valves used in the oil&gas industry often experience erosion due to the impingements of solids dragged by the extracted fluid. The considerable economic damages related to this phenomenon enhance the interest in techniques for predicting the location of the erosion hotsposts and the loss of material from the valve under aggressive flow conditions. CFD codes are often equipped with utilities for erosion estimation which rely on an Eulerian-Lagrangian model for two-phase flow computation and an erosion correlation for estimating the material removal caused by each particle-wall impingement. However, both the Eulerian-Lagrangian model and the erosion correlations include a number of uncertain submodels and parameterizations of empirical na...