The transport of solids in multiphase flows is common practice in energy industries due to the unavoidable extraction of solids from oil and gas bearing reservoirs. The persistent collision of solids to the pipeline can lead to erosion, i.e., the removal of internal surface of the pipeline. Reliable estimates of erosion rates are essential for designing and safely operating pipelines that transport solids. Prediction of erosion rates in multiphase flow is a complex problem due to the lack of accurate models for predicting particle movements in the flow and their impact velocities to the wall. The erosion-rate calculations also depend on the accuracy of the flow regime predictions in the pipeline. The comparisons of existing model predictions to experimental data revealed that the predictions might differ by several orders of magnitude for some operating conditions. The goal of this paper is to introduce a computational framework that estimates the model-prediction uncertainty of erosion-rate models. The inputs are a model predicting erosion rates and a database containing erosion-rate measurements at various operating conditions. The framework utilizes a non-parametric regression analysis, Gaussian Process Modeling (GPM), for estimating the model-prediction uncertainty. We compare two approaches for clustering the data prior to training GPMs: (1) a flow regime based clustering, and (2) a new clustering approach introduced in this paper. The results reveal that the new data clustering approach significantly shrinks the confidence intervals of the uncertainty estimates.