Uncertainty propagation is used to quantify the uncertainty in model predictions in the presence of uncertain input variables. In this study, we analyze a steady-state point-model for two-phase gas-liquid flow. We present prediction intervals for holdup and pressure drop that are obtained from knowledge of the measurement error in the variables provided to the model. The analysis also uncovers which variables the predictions are most sensitive to. Sensitivity indices and prediction intervals are calculated by two different methods, Monte Carlo and polynomial chaos. The methods give similar prediction intervals, and they agree that the predictions are most sensitive to the pipe diameter and the liquid viscosity. However, the Monte Carlo simulations require fewer model evaluations and less computational time. The model predictions are also compared to experiments while accounting for uncertainty, and the holdup predictions are accurate, but there is bias in the pressure drop estimates.Algorithms 2020, 13, 53 2 of 22 in Holm et al. [3,4], where they demonstrated how uncertainty propagation may be used for flow assurance on the Shtokman gas and condensate field. Their analysis included pressure drop and liquid holdup predictions using a one-dimensional model in the software OLGA. They determined probability distributions for a selection of input variables and closure laws and they propagated these uncertainties through the multiflow model using a Monte Carlo method. The result is the 10th, 50th, and 90th output percentiles for pressure drop and liquid holdup predictions and measures of sensitivity to the input variables. Hoyer et al. [5] used Monte Carlo simulations with OLGA in order to identify influential variables and closure laws in several groups of data with different flow conditions. They are only able to construct satisfactory probability distributions for each closure law when using a group of data where the closure law is significant.Klavetter et al.[6] modeled liquid holdup and pressure drop in two-phase pipe flow using the TUFFP Unified Model for two-phase flow. They assumed an uncertainty range for each input variable and compared perturbation, Taylor series approximations and Monte Carlo methods for uncertainty propagation. They concluded that Taylor series approximations overestimate the output uncertainty while the other methods perform well. Keinath et al. [7] also demonstrated the importance of selecting an appropriate framework when handling uncertainty in multiphase modeling and highlighted the value of quantitative information about the input uncertainty distributions for decision making. Just recently, in Liu et al. [8], a Gaussian process and principal component analysis were applied to a complex two-phase flow model in order to explore the uncertainty and reduce the complexity of the model. Picchi and Poesio [9] considered a one-dimensional model for two-phase pipe flow. Known distributions for input variables are propagated through the model using Monte Carlo methods to obtain first-order a...