This article presents a one-dimensional numerical model for vertical upward multiphase flow dynamics in a pipeline. A quasi-steady-state condition is used for the gas phase as well as liquid and gas momentum equations. A second-order polynomial for homogeneous flows and a sixth-order polynomial for separated flows are derived to determine the two-phase flow dynamics, assuming that the gas flow mass is conserved. The polynomials are formulated based on the homogenous and separate flows' momentum equation and the homogeneous flows' rise velocity equation and their zeros are the flow actual liquid holdup. The modeling polynomial approach enables the study of the polynomial liquid holdup zeros existence and uniqueness and as a result the design of a stable numerical model in terms of its outputs. The one-dimensional solution of the flow for the case of slug and bubble flow is proven to exist and to be unique when the ratio of the pipe node length to the time step is inferior to a specific limit. For the annular flow case, constraints on the inlet gas superficial velocity and liquid to gas density ratio show that the existence is ensured while the uniqueness may be violated. Simulations of inlet pressure under transient condition are provided to illustrate the numerical model predictions. The model steady-state results are validated against experimental measurements and previously developed and validated multiphase flow mechanistic model.