Most hydrocarbons production systems inevitably operate under multiphase flow conditions. While there is an agreement that Navier-Stokes equations govern variables for single-phase flow; there is no such consensus for the multiphase flow yet. Different numerical methods with dissimilar concepts are being conveniently used to simulate multiphase flow systems. Some of these methods do not respect the balance while others damp down strong gradients. The degree of complexity of these models makes the solution practically not reachable by numerical computations despite the fact that many rigorous and systematic studies have been undertaken so far. The essential difficulty is to describe the turbulent interfacial geometry between the multiple phases and take into account steep gradients of the variables across the interface in order to determine the mass, momentum and energy transfers. Different numerical techniques have been developed to simulate the gas-liquid simultaneous flow utilizing the CFD (Computational Fluid Dynamics) discipline. For example, the Volume of Fluid (VOF) model, the Eulerian-Langrangian models and the Eulerian-Eulerian models and combinations between these models. This paper presents the outcomes of numerical investigation carried out to probe the effect of a horizontal bend on the behavioral phenomenon of incompressible air-water simultaneous flow. A 3D CFD code has been developed based on NASA VOF code, which was designed for a different application. Major modifications were implemented on the original program to develop a fit-for-purpose one. The results have been qualified with the experimental data available from a different part of the same project and satisfactory agreements were obtained.