With the development of computer capability and advancement of modelling software, the technology of computational fluid dynamics (CFD) has been gradually recognized and utilized in engineering fields to help seeking practical industrial solutions of gas-liquid two-phase flow. Among various numerical models for two-phase flow, the two-fluid model which treats two phases separately by using two sets of transport equations can be considered as the most accurate one. However, the numerical descriptions of interfacial actions bridging these two phases become to critical aspects to transport equations. The interfacial process descriptions named closure terms generally include three aspects, referred as turbulence model, interfacial driving force and interfacial area concentration (interfacial area per unit). In current study, the capability of Two-fluid model and its closure terms has been investigated to simulate vertical and horizontal bubbly flow conditions.